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Abstract 



The production and decay of supersymmetric particles is presented in this thesis. The search 
for light mixing top squarks and neutralinos/charginos will be a major task at the upgraded 
Tevatron and the LHC as well as at a future e + e" linear collider. The dependence of the 
hadro-production cross section for weakly and strongly interacting particles on the renormal- 
ization and factorization scales is found to be weak in next-to-leading order supersymmetric 
QCD. This yields an improvement of derived mass bounds or the measurement of the masses, 
respectively, of neutralinos/charginos and stops at the Tevatron and at the LHC. Moreover, the 
next-to-leading order corrections increase the predicted neutralino/chargino cross section by 
+20% to +40%, nearly independent of the mass of the particles. The consistent treatment as 
well as the phenomenological implications of scalar top mixing are presented. The corrections 
to strong and weak coupling induced decays of stops, gluinos, and heavy neutralinos [including 
mixing stop particles] are strongly dependent on the parameters chosen. The decay widths are 
defined in a renormalization scheme for the mixing angle, which maintains the symmetry between 
the two stop states to all orders. The correction to the stop production cross sections, depending 
on the fraction of incoming quarks and gluons, varies between -10% to +40% for an increasing 
fraction of incoming gluons. The dependence of the cross section on all parameters, except for 
the masses of the produced particles, is in contrast to the light-flavor squark case negligible. The 
calculation of the stop production cross section was applicable to the Tevatron search for particles 
which could be responsible for the HERA anomaly. 



ZUSAMMENFASSUNG 

In dieser Arbeit wird die Produktion und der Zerfall supersymmetrischer Teilchen betrachtet. Die 
Suche nach leichten mischenden top-Squarks und Neutralinos/Charginos ist eine der wichtigsten 
Aufgaben am Tevatron und LHC ebenso wie an einem zukiinftigen e + e~-Linearbeschleuniger. 
Die Abhangigkeit der Wirkungsquerschnitte fur die Produktion stark und schwach koppelnder 
Teilchen von der Renormierungs- und Faktorisierungsskala ist in nachstfuhrender Ordnung re- 
duziert. Dies erlaubt verbesserte Massenschranken oder eine verbesserte Massenbestimmung fiir 
Neutralinos/Charginos und Stops am Tevatron und am LHC. Daruber hinaus vergroBern die Kor- 
rekturen den vorhergesagten Wirkungsquerschnitt fiir Neutralinos/Charginos um +20% bis +40%, 
nahezu unabhangig von der Masse der produzierten Teilchen. Weiterhin wird eine konsistente 
Behandlung der Stop-Mischung vorgestellt und deren phanomenologische Konsequenzen unter- 
sucht. Die Korrekturen zu starken und schwachen Zerfallen von top-Squarks, Gluinos und schw- 
eren Neutralinos — insofern sie top-Squarks enthalten — hangen significant von den gewahlten 
Parametern ab. Die Zerfallsbreiten enthalten die Definition des Stop-Mischungswinkels, welche 
die Symmetric zwischen den beiden Stop-Zustanden in beliebiger Ordnung Storungstheorie 
erhalt. Die Korrekturen zu den Wirkungsquerschnitten fiir Stop-Produktion variieren zwischen - 
10% und +40% und wachsen mit dem Anteil einlaufender Gluonen gegeniiber Quarks. Im Gegen- 
satz zur Produktion massenentarteter Squarks ist die Abhangigkeit von Parametern iiber die Stop- 
Masse hinaus vernachlassigbar. Die Berechnung des Wirkungsquerschitts fiir die Stop-Produktion 
konnte am Tevatron auf die Suche nach Teilchen, die fiir die HERA-Anomalie verantwortlich sein 
konnten, angewandt werden. 



Vernunft und Wissenschaft gehen oft verschiedene Wege 



Paul K. Feyerabend, Wider den Methodenzwang 
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Introduction 



A fundamental element of particle physics are symmetry principles. The electroweak as well as 
the strong interaction, combined to the Standard Model, are based on the gauge symmetry group 
SU(3)xSU(2)xU(l). The extension [[Ij] of this concept to a theory incorporating global or local su- 
persymmetry is a well-motivated step for several reasons: 

The Standard Model has been well-established by the discovery of the gluon and the weak gauge 
bosons, and by precision measurements at LEP and at the Tevatron, as well as at HERA. Currently, 
there is no experimental compulsion to modify the Standard Model at energy scales accessible to these 
colliders, provided the predicted Higgs boson will be found at LEP or at a future hadron or electron 
collider. However, a set of conceptual problems cannot be solved in the Standard Model framework: The 
mass of the only fundamental scalar particle, the Higgs boson, is not stable under quantum fluctuations, 
i.e. loop contributions to the Higgs mass term become large at high scales and have to be absorbed 
into the counter terms for the physical Higgs mass. This hierarchy problem leads to fine tuning of the 
parameters in the Higgs potential, to avoid the breakdown of perturbative weak symmetry breaking. 

Possible grand unification scenarios are based on a gauge group at some high unification scale, which 
contains the different Standard Model gauge groups. Simple unification groups are the SU(5) [Q] or 
SO(10) [|], the latter favored in scenarios with massive neutrinos. Non-minimal scenarios may yield 
intermediate symmetries and threshold effects, but as long as they include a simple unifying gauge group, 
the three running Standard Model couplings have to meet in one point at the unification scale. The 
requirement of one unification point and additional bounds from the non-observation of the proton decay 
lead to difficulties in the Standard Model, when it is embedded into a grand desert scenario, and most 
likely restrict the validity of the Standard Model to scales around the weak scale. 

In supersymmetric extensions of the Standard Model the masses of scalar particles remain stable even 
for very large scales, as required by grand unification scenarios. Quantum fluctuations due to fermions 
and bosons cancel each other; the leading singularities also vanish in softly broken supersymmetric 
theories. The hierarchy problem does therefore not occur in the extended supersymmetric Higgs sector. 
Including an intermediate supersymmetry breaking scale, the minimal supersymmetric extension of the 
Standard Model may be valid up to a grand unification scale without any fine tuning, being compatible 
with grand desert unification scenarios. Given the strong and the Fermi coupling constant at low scales, 
it predicts the weak mixing angle in very good agreement with the measured value M\. For a large top 
quark mass the renormalization group evolution can drive the electroweak symmetry breaking at low 
scales. The minimal supersymmetric Higgs sector consists of two doublets, in order to give masses to up 
and down type quarks while preserving supersymmeuy and gauge invariance. Hence, after breaking the 
weak gauge symmetry, five physical Higgs bosons occur. The non-diagonal CP even current eigenstates 
yield a light scalar Higgs boson with a strong theoretical upper bound on its mass. In some regimes of the 
supersymmetric parameter space this particle is accessible to LEP2, and the dependence of the theoretical 
mass bound on low-energy supersymmetry parameters can be used to constrain the fundamental mixing 
parameter tan (5. 
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In supersymmetric R parity conserving models the lightest supersymmetric particle is stable. This 
LSP, which in many scenarios turns out to be the lightest neutralino, is a possible candidate for cosmo- 
logical cold dark matter. 

In analogy to the gauge symmetries one may extend the global to a local supersymmetry. This 
invariance gives rise to higher spin states in the Lagrangean: a massless spin-2 graviton field and its 
spin-3/2 gravitino partner appear [0]. The general Einstein-gravitation is implemented into a theory of 
the strong and electroweak interaction. The so-obtained Kahler potential can in simple cases be derived 
by superstring compactification |^|. 

The breaking of exact supersymmetry is reflected in the observed mass difference between the Stan- 
dard Model particles and their partners. Due to the current mass limits, this mass difference is, in case 
of strongly interacting particles, much larger than the typical mass scale of the Standard Model particles. 
Assuming no mixing for light-flavor squarks, there are stringent mass limits on the squarks and gluinos 
from the direct search at the Tevatron [|7], H]. Due to large Yukawa couplings, the partners of the third gen- 
eration Standard Model particles may mix. Since more parameters of the supersymmetric Lagrangean 
enter through the non-diagonal mass matrices and the couplings, the mass limits for these third genera- 
tion sfermions are weakened. Moreover, all supersymmetric partners of the electroweak gauge bosons 
and the extended Higgs boson degrees of freedom mix. The search for these weakly interacting particles, 
neutralinos and charginos, at hadron colliders [Q] has not reached its limitations and will complete the 
limits obtained from the search at LEP2 The search for strongly interacting and also for light weakly 
interacting supersymmetric particles is one major task for the upgraded Tevatron and the LHC. The in- 
vestigation of mixing effects in the strong and weak coupling sector requires precision measurements at 
hadron as well as at lepton colliders. 

The reconstruction of supersymmetric particles from detector data is difficult in R parity conserving 
theories, since two LSPs leave the detector unobserved. Moreover, hadron colliders do not have an in- 
coming partonic state with well-defined kinematics, but the partonic cross sections have to be convoluted 
with parton density functions. The derivation of mass bounds or the mass determination, respectively, 
has to be performed by measuring the total hadronic cross section, if rather specific final state cascades 
cannot be used to determine the mass. Especially for strongly interacting final state particles, the cross 
sections depend on the factorization and renormalization scales through the parton densities and the run- 
ning QCD coupling. The scale dependences lead to considerable uncertainties in the determination of 
mass bounds. The next-to-leading order cross sections will improve the mass bounds not only by their 
accuracy but also by their size. These hadronic cross sections for mixing supersymmetric particles at the 
upgraded Tevatron as well as at the LHC will be given in this thesis. 



Similarly to light-flavor squarks and gluinos, the search for top squarks with a non-zero mixing angle 
will lead to stringent mass bounds, which are essentially independent of the mixing parameters and the 
masses of other supersymmetric particles. However, it will most likely be impossible to measure the 
mixing angle at hadron colliders directly, since the cross sections for the production of a mixed stop pair 
are strongly suppressed. The analysis of mixing effects in the stop sector will be completed by the direct 
measurement of the mixing angle in e + e~ collisions [10]. 

Regarding certain decay channels the direct search for gauginos and higgsinos at hadron colliders 
resembles the search for weak gauge bosons. Although in most supergravity inspired scenarios not all 
gauginos and higgsinos are light enough to be found at the upgraded Tevatron, the search for light neutral 
and charged gauginos is promising and could improve the LEP2 results at the upgraded Tevatron and at 
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the LHC. Even if the leading order cross sections are independent of the QCD coupling, they depend 
on the factorization scale through the partem densities. The next-to-leading order predictions will again 
considerably improve the bounds derived for masses and couplings. 

However, all search strategies for supersymmetric particles depend on cascade decays leading to 
leptons, jets, and LSPs in the final state, the latter provided R parity is conserved. As long as the masses 
of the particles under consideration are not known, the analysis of these multiple decay channels does 
not give strong limits, e.g. on mixing parameters involved. But for a sufficiently large sample of events 
including supersymmetric particles, the whole variety of possible decays and couplings will help to 
determine the mass and mixing parameters of the supersymmetric extension of the Standard Model. The 
measurement of low energy parameters can then be used to search for universal parameters, predicted by 
grand unification or supergravity inspired scenarios. 



Outline of the Thesis 

Since the supersymmetric observables presented in the following analyses can, from a phenomenological 
point of view, be treated independently, technical features are covered with their first appearance. 

The general physics background is described in the first chapter. A short introduction into supersym- 
metric extensions of the Standard Model is complemented by the discussion of special aspects concerning 



mixing particles; the next-to-leading order treatment of the mixing angle [11] in the CP conserving stop 
sector is presented, and the regularization prescriptions used for supersymmetric gauge theories are sum- 
marized. The supersymmetric Feynman rules and a complete set of formulae considered useful for the 
detailed understanding of the calculations are given in the appendices. 

The production cross sections for neutralinos and charginos at hadron colliders are treated in chap- 
ter ^. They include the virtual and real next-to-leading order corrections, the latter calculated using the 
dipole subtraction method. The treatment of on-shell singularities is described in detail. The possible 
improvement of the current analysis by using the next-to-leading order cross section is pointed out. 

In chapter || the decay widths including mixing stop particles in next-to-leading order supersymmetric 
QCD [ |i"l| ] are given. They include weak and strong coupling stop decays as well as gluino and heavy 
neutralino decays to a light stop. The treatment of the mixing angle follows the theoretical description in 
chapter [T[ The complete analytical results for the next-to-leading order stop decay width is presented in 
the appendix. 

The study of scalar top quarks is continued in chapter |], where the production cross section at hadron 



colliders is given for both of the mass eigenstates in next-to-leading order supersymmetric QCD [12]. 
One crucial point is the influence of the mixing angle and supersymmetric parameters, which are present 
in the virtual collections, on the experimental analysis and on the mass bounds. The real gluon emission 
is calculated using the cut-off method. 

For a light stop the production cross section at hadron colliders can be adapted to R parity violating 
scenarios [|i"4]]. The resonance cross section for the production of R parity violating squarks in ep colli- 



sions is calculated in next-to-leading order [13], and the search results for these particles at HERA and 
at the Tevatron are combined. The influence of other search strategies for R parity violating squarks is 
reviewed. 

The analytical calculations have been performed using the symbolic manipulation program 
FORM [15], for the numerical integration routine VEGAS [jl^] was chosen, and the parton cross sections 



have been calculated using the CTEQ4 [IT7p parton densities in leading and next-to-leading order. 



1. SUPERSYMMETRY 



1.1. Supersymmetric Extensions of the Standard Model 

Global supersymmetry is a possible extension of the set of symmetries appearing in flat space-time gauge 
theories. The most general extension of a Poincare invariant theory would be an iV-extended super- 
Poincare Algebra containing central charges Jl8|, The supersymmetry generators Q % [i = 1,...,N] 
and their complex conjugate Q 1 transform fermionic into bosonic fields and vice versa, therefore obey- 
ing an anticommutation relation. These anticommutators lead to a Z2 graded Lie algebra, containing the 
supersymmetry as well as the Poincare group generators and circumventing the No-Go theorem [ ]. 



The dimension of the extension N determines the maximum spin present in the particle spectrum of 
the theory. Renormalizability requires a maximum spin of one for global supersymmetry, which is 
equivalent to N < 4. Including the graviton results a maximum spin two for local supersymmetry, 
supergravity, and renders N < 8. For N = 1 this super-Poincare algebra becomes particularly simple, 
since the central charges vanish and the generators Q, Q anticommute with themselves. Extended super- 
symmetric theories have some remarkable features: for N = 2 the particle spectrum can be calculated 
non-perturbatively [fEOQ, N = 4 leads to a completely finite theory, and N > 5 contains gravitation. 
However, the observed low energy particle spectrum and CP violation are only compatible with (N=l) 
global supersymmetry. We will make use of the incorporation of global supersymmetry into local super- 
gravity only by assuming certain characteristics of the mass spectrum at high scales, where unification is 
required. 

Since supersymmetric theories by definition contain scalar particles not only in the Higgs sector, the 
behavior of scalar masses is of importance: In the Standard Model the scalar Higgs boson mass suffers 
from UV divergent radiative corrections, proportional to g 2 A. 2 where g is a gauge coupling and A is an 
UV cut-off parameter. This cut-off parameter could be fixed by some scale where new physics appears. 
Assuming the Standard Model not being an effective theory for mass scales around the weak gauge boson 
mass, e.g. leads to a physical scalar mass of the order of the weak scale and higher order loop contribu- 
tions of the order of the cut-off, which could be the Planck scale. These corrections have to be absorbed, 
using fine-tuning of mass and coupling counter terms in the Lagrangean. The large corrections in the 
Standard Model originate from gauge boson and top quark loops. In supersymmetric extensions addi- 
tional collections arising from the supersymmetric partners enter with a minus sign and weaken the UV 
degree of divergence to a logarithmic behavior [5m /m oc log A 2 ]. For broken supersymmetry another 
term proportional to the mass difference between the Standard Model loop particles and their supersym- 
metric partners arises. Assuming e.g. a grand desert SU(5) scenario^] the natural shift of the scalar masses 



'Any Lie group containing the Poincare group and a compact inner symmetry group factorizes, i.e. the generators of the 
Poincare group and the inner symmetry group commute with each other. The extended Lie algebra becomes trivial. 

2 Though not all matter fields can be unified in one SU(5) multiplet. A more generic GUT model would be supersymmetric 
SO(10), directly broken to the Standard Model gauge group. In contrast to SU(5), SO(10) unification with non-zero neutrino 
masses may lead to the observed baryon asymmetry [M], The numerical analyses e.g. of gauge coupling unification in grand 
desert SU(5) and SO(10) scenarios are similar. 



H [GeV] 

Figure 1.1: The renormalization group evolution of the three Standard Model gauge couplings, assuming 
a simple SU(5) GUT gauge group and a grand desert [22 ]. 



between the weak and the unification scale is limited to less than one order of magnitude [22]. 



Assuming a supersymmetric extension of the Standard Model, the evolution of the gauge couplings 
can be evaluated, based on different scenarios. Embedding the Standard Model into a simple GUT gauge 
group does neither fix the gauge group nor possible intermediate scenarios, i.e. SU(5) unification with a 
grand desert is only one possibility. The evolution of the gauge couplings depends on threshold effects 
and masses in intermediate unification models. However, it can be shown, that, using supersymmetry, 
the gauge couplings unify up to a certain accuracy, Fig. |1.1| , The unification of the three Standard Model 
gauge couplings determines one of the three parameters involved [a, = sin 2 9 W , a s ] theoretically. 
This prediction has to be compared to the measured value e.g. for the weak mixing angle in the MS 
scheme. In contrast to the Standard Model, which for reasons described above will hardly be valid up to 
a large unification scale, the predicted value for the minimal supersymmetric extension of the Standard 
Model s 2 w {m z ) = 0.2334 ± 0.005 agrees very well with the measured value of 0.2316 ± 0.0003 [Q]. 
Any specific GUT scenario fixes the renormalization group evolution of all the masses and couplings. 
The determination of a s {m z ) from a and again reflects the improvement of the Standard Model one- 
scale GUT, which gives a s (m z ) = 0.073 ± 0.002, to the supersymmetric one-scale GUT, which yields 
a s {m z ) = 0.129 ± 0.010. However, the measured value of a s (m z ) = 0.118 ± 0.004 indicates, that the 
minimal supersymmetric GUT prefers a slightly larger value; threshold effects may be responsible for 
the difference. Very light gauginos and very heavy squarks might even for the SU(5) GUT model lead to 



the measured value of a s [22]. 



1.2. The Minimal Supersymmetric Standard Model 

The Minimal Supersymmetric Standard Model [MSSM] adds a minimal set of supersymmetric partner 
fields to the Standard Model [SM]. These fields contain scalar partners - sleptons and squarks - of all 
chiral eigenstates of the Dirac fermions, incorporated into chiral supermultiplets. Absorbing Standard 
Model gauge fields into vector supermultiplets leaves Majorana fermion^] partners of the neutral U(l), 

3 Majorana fermions are denned as their own anti-particles, i.e. the Majorana spinor is constructed by combining two Weyl 
spinors. Some arbitrariness may arise from the different treatment of electric and color charge, which leads to Majorana 



6 



CHAPTER 1. SUPERS YMMETRY 



SU(2), SU(3) gauge fields and Dirac fermion partners of the charged SU(2) gauge fields, called gauginos. 
The SU(3) ghost fields are defined by re-writing the Fadeev-Popov determinant; therefore they do not 
receive supersymmetric partners, which enter by requiring the original Lagrangean being invariant under 
a global supersymmetry transformation. 

For reasons described later in this chapter the supersymmetric scalar potential cannot include conju- 
gate fields. Hence, at least two different complex Higgs doublets have to be introduced to give masses 
to up and down type quarks. They form five physical Higgs particles after breaking SU(2)xU(l) in- 
variance. This extension of the Higgs sector is not generically supersymmetric. However these scalar 
Higgs degrees of freedom have to develop partner fields. This yields neutral and charged Majorana/Dirac 
fermions with the same quantum numbers as the SU(2) gauginos. 

The supersymmetric Lagrangean in superspace can be constructed by extending the integration over 
the Lagrange density from space-time to a superspace integration, i.e. by adding two Grassmann di- 
mensions {8, 8}. All superfields can be written as a finite power series in these Grassmann variables, 
containing the component fields in the coefficients. Only two elements enter the supersymmetric La- 
grangean: (i) the so-called F term of a chiral supermultiplet, denoted as in the expansion of the 
superfield in the Grassmann variable 8; (ii) the D term of a vector multiplet Vq2§2 [jl|]- The kinetic 
real vector supermultiplet is defined as the product of a chiral supermultiplet and its conjugate Qj&j, 
its D terms contain the F components of the chiral multiplets F*Fj, which is absorbed into the scalar 
potential. 

The most general ansatz for a superpotential formed by chiral supermultiplets relies on the fact, that 
the product of two chiral supermultiplet is again chiral: 

W({$}) = mij^j + A, / 7, < i> ( «i > ; <i>/, (1.1) 

Higher orders in the polynomial would give mass dimensions bigger than four and therefore spoil renor- 
malizability. The superpotential occurs in the Lagrangean as (W + W). The scalar potential in the 
component-field Lagrangeanf] contains, after integration over the Grassmann variables, the non- Yukawa 
terms arising from the superpotential W; it is defined as 



dW{A) dW{A ) 
dAj j + d~Aj 



V = - ( F*Fj + —^ Fj + ~^F* ) (1.2) 



The Euler-Lagrange equations yield F* = — dW(A)/dAj where Aj are the sfermion fields in the su- 
permultiplet. This fixes the most general scalar potential including chiral supermultiplets in the matter 
sector of the Lagrangean: 

V = J2\Fjf (1-3) 



Including the gauge sector for a non-abelian gauge group leads to vector multiplets containing the 
gauge fields and their partners. The scalar potential will also contain the D auxiliary component field 

neutralinos and gluinos, but Dirac charginos. But it is just a name for the particles. Dirac charginos also yield fermion number 
violating vertices. 

4 A product of a chiral superfield and a conjugate is not a chiral but a vector multiplet, as the kinetic superfield. The su- 
perpotential therefore does not contain conjugate superfields and neither does the scalar potential contain conjugate component 
Higgs fields. 
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terms of the gauge multiple^] 

v = E i^-i 2 + \ E (^ a ) 2 = E i^-i 2 + y E (5*™) 2 a-4) 

3 a 3 a 

The D term is written for a general non-abelian SU(N) gauge group. S are the scalar fields transforming 
under the fundamental representation of the corresponding gauge group, and T a the generators of the 
underlying gauge group. 



1.2.1. R Parity 



The most general superpotential as given in eq.(l.l) contains trilinear couplings of chiral matter super 



multiplets, like the Higgs, the quark, and the lepton supermultiplet. Couplings between the different 
Higgs fields or between the Higgs field and corresponding lepton or quark supermultiplets are needed 
to construct the two doublet Higgs sector in the scalar potential. Although these vertices conserve the 
over-all fermion number, they may violate the baryon and lepton number and would lead to the same 
effects as leptoquarks, e.g. proton decay [|||]. In extensions of the Standard Model these operators are 
forbidden by gauge invariance, as long as their dimension is less than six. The MSSM either needs to 
suppress the different couplings or remove the whole set by applying a new Z2 symmetry which changes 
the sign of the Grassmann variables in the Lagrangean. The corresponding conserved charge is defined 
as 

R = (_1)3B+L+2S (1 5) 

where B is the baryon number, L the lepton number, and S the spin of the particle. This number is chosen 
to give (+) for Standard Model particles and (— ) for supersymmetric partners. The Higgs particles in 
the two doublet model are all described by R = +1. Accounting for R symmetry in the supersymmetric 
Lagrangean removes trilinear chiral supermultiplet vertices containing no Higgs superfield. The general 
superpotential eq.(|l.l[) can be separated into an R parity conserving and an R parity violating part, which 
read for one generation of quarks and leptons 

w = w R + 

W R = \ E EH{L j + \ D DH{Q j + X u UH J 2 Q j - ^h{h{ 

= -\VVE + X'L j Q j D + -\"UDD (1.6) 

The contraction of two indices is defined by the antisymmetric (2 x 2) matrix e^ ; L, Q are electron and 
quark doublet superfields, E, D, U are the singlet superfields for the electron, d and u type quark; A ' ,u 
are the Yukawa coupling matrices, and fi is the Higgs mass parameter, which also defines the higgsino 
mass [see appendix The Yukawa couplings A, A' violate lepton number, A" violates baryon number. 

The combination (A' • A") leads to proton decay via an s channel d type leptoquark and therefore has 
to be strongly suppressed. The conservation of R is therefore a sufficient condition for the stability of the 
proton. However, this symmetry has been introduced ad hoc for weak scale supersymmetry as a less rigid 
substitute for the conservation of some combination of B and L. The exact vanishing of (A' • A") is not 



~ The supermultiplet constructed from the gauge vector multiplet and including the field strength component field is chiral. 
But its F term contains the D component fields of the vector gauge multiplet. 
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a necessary condition for a stable proton; i.e. if is not removed by hand by demanding R symmetry 
in the supersymmetry Lagrangean, then many different constraints can be imposed on combinations of 
couplings and masses in the R parity violating sector. The limits from direct production at HERA as well 
as from rare decays generically determine A • \/m 2 , dependent on the flavor of the squark considered. 
The same holds for atomic parity violation. The bounds from neutral meson mixing influence A • A/m, 
and the direct searches at hadron colliders and LEP are only sensitive to the mass, except for the analysis 



of specific decay channels [24, 25] 



Phenomenologically, exact R parity conservation leads to the existence of a stable lightest super- 
symmetric particle (LSP), and allows for the production of supersymmetric particles only in pairs. For 
cosmological reasons this LSP has to be charge and color neutral, which restricts the choices in the 
MSSM framework to the lightest neutralino or the sneutrino. In GUT models exact R parity conserva- 
tion is not necessary to obtain low-energy R parity conservation. For broken R parity the unstable 'LSP' 
could therefore be long-living and charged, allowing for charginos, sleptons and even stops, as long as 
the lifetime is small enough to circumvent the cosmological constraints. 

1.2.2. Soft Breaking 

If supersymmetry would be exact, the squarks and sleptons were mass degenerate with the Standard 
Model particles. Since the gauge couplings have to respect supersymmetry in order to cancel the 
quadratic divergences, breaking supersymmetry means enforcing a mass difference between Standard 
Model particles and their supersymmetric partners. The mechanism of introducing mass terms by soft 



breaking [26] at a given scale has to respect gauge symmetry, weak-scale R parity, stability of scalar 
masses, and experimental bounds e.g. on FCNC. Soft breaking terms can be added to the superpotential 
eq.( |l.l[ ) at any given scale. They exhibit the generic form 

6' 



Asoft = - [n4) i:j C*Cj - - (m 1/2 ) . XjXj + h.c. 



-AijkdCjCk + B^HiH 2 + h.c. 



(1.7) 



The component fields involved are generic scalars C, Majorana fermions A and the Higgs fields Hi, H 2 , 
which are again contracted using gy. The possible set of parameters consists of: 

- Scalar mass matrices (m^j = 1, ...n] for squarks and sleptons with n generations. The 
diagonal masses can be chosen real, since (W + W) enters the Lagrangean. 

- Three real gaugino masses (mi/2)j [j = 1, 2, 3]. 

- 27 complex trilinear couplings Aijf.[i, j, k = 1, 2, 3] which conserve the R charge. 

- Two masses for the Higgs scalars and a complex Higgs mass parameter \iBH\ H 3 2 . 

Evolving soft breaking mass terms by means of the renormalization group equations can lead to breaking 
of the U(l)xSU(2) symmetry by driving one mass squared negative. This generalization of the Coleman- 



Weinberg mechanism p7| ] links the large top Yukawa coupling to electroweak symmetry breaking. 

1.2.3. Supersymmetric QCD 
Particle Content 



The search for directly produced supersymmetric particles at hadron colliders is dominated by strongly 
interacting final states. In these production processes the quantum corrections in next-to-leading order 
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are expected to be significant. Moreover, the corrections to the production of weakly interacting particles 
at hadron colliders are dominated by strong coupling effects. Although the parton picture and thereby the 
incoming state is not affected by the heavy supersymmetric partners of quarks and gluons, a consistent 
description of virtual particle effects requires the inclusion of these particles. 

The supersymmetric extension of the QCD part of the Standard Model is straightforward, since the 
SU(3) invariance is unbroken. One chiral mass superfield Q contains the left handed quark doublets 
(u L , d L ) and their squark partners (u L , d L ). Two more superfields U, D connect the quark singlet fields 
(u c R , d c R ) to their partners (u* R , d* R ). The SU(3)c x SU(2)l x U(l)y quantum numbers for quarks and 
squarks are identical. Whereas the q L is a SU(3) triplet, the q R is an anti-triplet and couples with (— T a ) 



to the quark and gluino, as can be seen in Fig. A.2. The gluon vector superfield mirrors the gluons to 



gluinos (g), which are real Majorana fermions and therefore carry two degrees of freedom^. The number 
of generations is not restricted by supersymmetry. The CKM matrix for the quarks will in the following 
be assumed to be the unity matrix. The same holds for the squark CKM matrix, which is not fixed by 
first principles to be either diagonal or equal to the quark matrix. 
The general mass matrix for up-type squarks is given by 



\ -m q (A q + //cot/3) m^+mq + |s^,m|cos(2/3) J 

For down type squarks cot (3 in the off-diagonal element has to be replaced by tan j3. The entries 
mQ,mu, A q are the soft breaking masses. In the diagonal elements the quark mass still appears, as 
in exact supersymmetry. The m z contributions arise from the different SU(2) quantum numbers of the 
scalar partners of left and right-handed quarks. For light-flavor squarks this matrix can be assumed being 
diagonal, since the chirality flip Yukawa interactions are suppressed. For the top flavor these off-diagonal 
elements cannot be disregarded. Taking into consideration bottom-tau unification the ratio of the Higgs 
vacuum expectation values tan j3 has to be either smaller than ~2.5 or larger than ~40. In the second 
case, a large value for tan compensates for the small bottom quark mass and yields a strongly mixing 
sbottom scenario. The results for the stop mixing may be generalized to the sbottom case. 

Neglecting additional mixing from a CKM like matrix, the chiral squark eigenstates are equal to the 
mass eigenstates for the light flavors. If we furthermore assume the soft breaking mass being dominant 
and invariant under SU(2), then the light-flavor mass matrix is proportional to the unity matrix, i.e. the 
masses of the ten light flavor squarks are equal. As long as only strong coupling processes are considered, 
we will have to deal with ten identical particles. This will not be the case for the scalar top sector as will 
be shown in section |1.2.4 



m 



+m 2 q +(\ - \ s w) m l cos(2/3) —m g (A g + ^cot 0) 



1.2.4. Mixing Stop Particles 
Diagonalization of Mass Matrices 



For scalar top quarks the off-diagonal elements of the squark mass matrix eq.(|1.8|) are large. Any real 
symmetric mass matrix of the form 



M\ L M\ r 
M 2 LR M\ R 



(1.9) 



The matching of the degrees of freedom is a subtlety in dimensional regularization, see section 



1.5 
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g 



Figure 1.2: Feynman diagrams for the stop self energy in NLO, including mixing in the second and third 
diagram, the top-gluino loop and the pure squark tadpole. 



can be diagonalized by a real orthogonal transformation, i.e. a uniquely denned real rotation matrix. The 
eigenvalues are 

o 1 



m. 



Tx{M 2 ) T [Tt 2 (M 2 ] - 4Det(7W 2 )] 



2MV2 



12 2 

The cosine of the mixing angle can be chosen positive — 7r/4 < 8 < 7r/4: 



cos(2<9) 



\Ml L -M 2 RR \ 



Tt 2 {M 2 ) -4Det(M' 



sin(20) 



Tr 2 (M 2 ) -4Det(>f 2 ) 



(1.10) 



(1.11) 



There is no flat limit from different to equal mass eigenvalues for this diagonalization procedure, since the 
diagonalized matrix would be proportional to the unity matrix and therefore commute with any rotational 
matrix. 



Stop Mixing 

In the scalar top sector the unrenormalized chiral eigenstates are t L o,t R Q. The chirality-flip Yukawa 



interactions give rise to off-diagonal elements in the mass matrix eq.(|1.8|) i.e. the bare mass eigenstates 
iio an d ho are obtained by a leading-order rotation, as described above. 



ho \ = ( cos6 sin § \ ( ho 
ho J V -sin 0o cos 6*o / V ho 



(1.12) 



The mass eigenvalues and the leading-order rotation angle 9q can be expressed by the elements of the 
mass matrix. However, SUSY-QCD corrections, involving the stop and gluino besides the usual particles 
of the Standard Model, modify the stop mass matrix and the stop fields. The Feynman diagrams are given 
in Fig. 1.2. As described in appendix ^, the coupling to a quark and a gluino as well as the coupling 



between four squarks can switch the chirality state and therefore contribute not only to the diagonal but 
also to the off-diagonal matrix elementsTI. This gives rise to the renormalization of the masses and of 

1/2 ~ 

the wave functions [£^o = Z^ tj]. Any leading-order observables concerning the mixing top squarks 



are linked by a re-rotation of ir/2, denoted by V\2, eq.( |A.2| ). In next- to-leading order this symmetry 
is broken by the mixing stop self energy. In order to restore this symmetry in any order perturbation 
theory^, we choose a real wave-function renormalization matrix Z 1 / 2 , which is defined to split into a 

7 It can be shown that a correction to the mass matrix renders the NLO mass matrix complex symmetric and not hermitian, 
as long as CP is conserved, i.e. only imaginary parts from the absorptive scalar integrals arise. 

8 Any observable containing only one kind of external stop particles can be transformed by exchanging the stop masses 
and adding (-) signs to sin(2#) and cos(26 l ). This prescription w ill be used for stop decay widths and for the hadronic 
production cross section in LO and NLO later and is defined in eq.(|A.2[), 
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real orthogonal matrix 71(89) and a diagonal matrix Z^^L, i.e. Z 1 / 2 = 71(89) Z X P 1 . The rotational part 



can be reinterpreted as a shift in the mixing angle QlT] , [28| ], given by 9 — 59 = 9: 

(1.13) 



f \ / 7~ 1/2 n 

*1 \ / ^diag.ll " 




This counterterm for the mixing angle allows the diagonalization of the real part of the inverse stop 
propagator matrix in any fixed-order perturbation theory. 



Re[^(p 2 )] =(Z 1/2 ) T b 2 l--M 2 + ReS(p 2 )] (Z^) 

= (4a g ) T [p 2 1 " n^r 1 {M 2 + ReS(p 2 )) 11(59) 
= (4Cg) Re^Jp 2 ) (4(g) (1.14) 



(4/ag) 



This holds as long as the real part of the unrenormalized stop self-energy matrix ReS(p 2 ) and thereby 
the whole next-to-leading order mass matrix is symmetric^ The mixing angle depends on the scale of 
the self energy matrix 

, f oxSf 2 A 2ReS 12 (p 2 ) 2ReS 12 (p 2 ) 2a 

V / mf — m~ + ReL22(P ) — ReLii^) to- — m~ 

We fix the renormalization constants by imposing the following two conditions on the renormalized stop 
propagator matrix: (i) the diagonal elements should approach the form l/D Ten jj(p 2 ) — > p 2 — m~ + 

im~ t Xi. for p 2 — > m| , with mj. denoting the pole masses; (ii) the renormalized (real) mixing angle 9 

is defined by requiring the real part of the off-diagonal elements A-en, 12 (p 2 ) and -D ren , 21 (p 2 ) to vanish. 
The three relevant counter terms for external scalar particles are 

K = ReEtf(m? ) i% = -ReS,,(m| ) ^(p 2 ) = - ^g^T d-17) 

Thus, for the fixed scale p 2 the real particles t\ and i 2 propagate independently of each other and do not 
oscillate. 

The so-obtained (running) mixing angle depends on the renormalization point Q, which we will indi- 
cate by writing 9(Q 2 ). The appropriate choice of Q depends on the characteristic scale of the observable 
that is analyzed. The real shift connecting two different values of the renormalization point is given by 
the renormalization group, leading to a finite shift at next-to-leading order SUSY-QCD 



9(Ql) - 9(Q 2 2 ) = CFasm l mt C ° 2 s(2g) Re[ B(Q 2 - m- g , m t )-B(Q 1 ; m~ g ,m t )] (1.18) 



9 The next-to-leading order SUSY-QCD correction to the stop mass matrix is 



Si2(p 2 ) = -2-kCfo.s ys A gA(m^) - s 4 gA(m {l ) + 8m§m t c 2 gB(p;m§, m t )J = £ 2 i(p 2 ) 

Sn(p 2 ) = -4nC F a s [(l +clg)A(m h ) + s 2 2§ A(m i2 ) - 2A(m- g ) - 2A(m t ) 

- 2(p 2 + m 2 1 )B(p;X,m il ) + 2{p 2 - m§ - m 2 + 2m g m t s. 2§ )B{p;m g ,m t )^ = P 1 " 2 1 E 22 (p 2 ) 

E li (p 2 ) = aE l3 (p 2 )/9p 2 (1.15) 
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200 400 600 800 1000 

Q [GeV] 

Figure 1.3: The dependence of 0(Q 2 ) on the renormalization scale Q. The input mass values are: 
m i/2 = 150 GeV, rriQ = 800 GeV, Aq = 200 GeV, p, > 0, for which the leading-order mixing angle 
is given by 1.24 rad. The minimum of the correction corresponds to the threshold Q = m g + mt in the 
scalar integral. 



This shift is independent of the regularization. In the limit of large scales the difference behaves as 
\og{Q 2 /Q' 2 ). A numerical example is presented in Fig L3. As a noteworthy consequence of the running- 



mixing-angle scheme, we mention that some LO symmetries of the Lagrangean are retained in the NLO 
observables. For instance, if for only one kind of external stop particle one chooses Q = m t ~, the results 



for the other stop particle can be derived by the simple operation V12, eq.(|A2j), which then also acts on 
the argument of the mixing angle. 

Considering virtual stop states with arbitrary p 2 , the off-diagonal elements of the propagator matrix 
can be absorbed into a redefinition of the mixing of the stop fields, described by an effective (complex) 
running mixing angle c g(p 2 ) = 9 — 59 c fi(p 2 ). This generalization amounts to a diagonalization of 
the complex symmetric stop propagator matrix D ren , including the full self-energy £(p 2 ), by a complex 
orthogonal matrix 7£(5# e fj)[^] exactly in analogy to eq.(1.14). The so-defined effective running mixing 
angle is given by 



9 eS (p 2 ) 



9n arctan 

u 2 



2E 12 (p 2 



m- 



m| + £22 (p 2 



£11 (p 2 



NLO 



9( P 2 ) + 



ImSi 2 (p 2 



ra- 



ti 



m~ 



(1.19) 



The complex argument of the trigonometric functions leads to hyperbolic functions. From this point 
of view the use of a diagonal Breit-Wigner propagator matrix is straightforward. For instance, in the 
toy process tg — ► tg all NLO stop-mixing contributions to the virtual stop exchange can be absorbed 
by introducing the effective mixing angle in the LO matrix elements. The argument of this effective 
mixing angle is given by the virtuality of the stop particles in the s channel. This procedure also applies 



Any real symmetric matrix can be diagonalized by a real orthogonal transformation O AO where O 



O. One 



generalization is the complex unitary diagonalization of a complex symmetric matrix U T AU with IP — U~ L , where the 
diagonal matrix is real and positive. Another one is the complex orthogonal diagonalization of a complex symmetric matrix 
O t AO, O -1 = O t where the diagonalized matrix is still complex. Note that a hermitian matrix can only be diagonalized by 
a unitary transformation U -1 AU. 
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to multi-scale processes like qq — ► tUg or e + e — > iii2, where the effective gi^/ji^ couplings 
become non-zero due to the different scales of the redefined stop fields. 

There exist other renormalization schemes for the stop mixing angle, either fixing the scale of the 
running mixing angle at some appropriate scale or absorbing certain diagrams e.g. contributing to the 



production process e + e~ — > tit 2 [|29p. Any of these schemes can be regarded as a prescription to mea- 
sure the mixing angle, either in the mixed production at e + e _ linear colliders or in decay modes or 
quantum corrections. The mixed production induced scheme however has the disadvantage of introduc- 
ing the Zt\t2 weak coupling constants into the QCD counter terms. The measured values of the mixing 



angle can be translated from one scheme into another by comparing the counter terms. In Fig. 1.3 the 
numerical effect of the finite renormalization can be seen to be small; the same holds for the different 
renormalization schemes, which are numerically almost equivalent. 

When fixing the counter term for the stop mixing angle 6, one can express the angle in terms of 



the parameters appearing in the mass matrix eq.(|L8|). The counter term 5(sin(20)) can be linked to the 
counter terms of these parameters: 

2m t (A t + jucot/?) 



b 2§ ~ 


2 


2 




- 

ti 


- m~ 




Ss 26 _ 


5(m~ 
ti 


— TO? ) 


Suit 
+■ — - 


S 28 


2 

to- 
ti 


2 

— m~ 
ti 


m t 



+ 5(A t + ncot f3) 
A t + //cot/3 

where 5x denotes the counter term of the parameter x. Since [i and appear in the scalar potential 
only in the weakly interacting sector, they will not be renormalized in next-to-leading order SUSY-QCD. 
However, 6 At can be calculated from the mass and mixing angle counter terms. This reflects the fact, 
that the system of observables used in the Feynman rules is non-minimal, i.e. the on-shell scheme for the 
masses and the running mixing angle determine the renormalization of the couplings ii^G and tifyA , 



where A t appears explicitly [30]. 



1.3. GUT inspired Mass Spectrum 

Next-to-leading order calculations in the framework of light-flavor SUSY-QCD ^ only incorporate a 
few free parameters: the Standard Model set and the gluino and the light-flavor squark mass. Including 
mixing stops and the mixing neutralinos/charginos the number of low-energy parameters becomes large. 
Hence, for a rough phenomenological analysis we will use a simplifying scenario, which could be a 



SUSY-GUT scenario, either supergravity [H] or gauge mediation [31] inspired. 



SUSY-GUT Scenario 

Inspired by the unification of the three Standard Model gauge couplings in supersymmetric GUT models 
we will assume a relation between these couplings and the gaugino masses. Independent of the actual 
form of the simple gauge group and the connected GUT scenario, the three Standard Model gauge groups 
are embedded into, and independent of intermediate scale particles and thresholds, we can assume gauge 
coupling unification. 

Mi(Q) _ M 2 (Q) _ M 3 (Q) _ rn 1/2 (M x ) 
a\{Q) a 2 (Q) a 3 (Q) a G u T (M x ) 
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where my 2 is the mass entry in the scalar potential, defined at the unification scale Mx- There the three 



gauge couplings unify to oqut — 1/26. For the masses at the weak scale this leads to J33| ] 



nig ~ M\{mz) — 0.4 mu 2 
mw - M i( m z) - 0.8 m 1/2 

m-g ~ M 3 (m z ) ~ 2.6 m l/2 (1.22) 

However, the gluino mass is strongly dependent on the scale which can lead to a difference of 30% 
between the pole mass m g and the running mass M 3 (M 3 ) For the derivation of these mass relations 
it is only necessary to assume a simple unification gauge group arising at a scale Mx ~ 2 • 10 16 GeV. 
The gaugino mass unification can be tested experimentally at the LHC [|2|] as well as at a future linear 



collider [10]. 



Mass Unification 

In a supergravity inspired MSSM [mSUGRA] the scalar masses and the trilinear couplings are assumed 
to be universal at the unification scale MxQ- In simple supergravity models they depend on the gravitino 
mass scale m 3 / 2 and on the cosmological constant [§]. The universal parameters at the unification scale 
Mx will be refered to as tuq and Aq. The parameter [iB occuring in the Higgs sector of the scalar 



potential [section |1.2.2| ] will be fixed by the choice of mi/2, mo, Ao, tan/3 and the Standard Model 
parameters, and by the requirement of electroweak symmetry breaking, up to its sign. The light-flavor 
squark masses can be expressed in terms of the universal scalar and gaugino masses, the other parameters 



only enter the off-diagonal elements of the mass matrix eq.(1.8) and can be neglected 



m\ h ~ ml + 6.3m? /2 + 0.35D m\ R ~ m\ + 5.8m 2 1/2 + 0.16-D (1.23) 

where D = m 2 z cos (2/3) < 0. For mSUGRA scenarios a general prediction for the light-flavor squark 
mass can be given [ ]33| ] 

mq > 0.85m g (1.24) 

Approximate Solution 

The stop masses can be expressed in terms of the top Yukawa coupling Y t = /i|/(47r). For small tan (5 
they approximately read 

M 2 f-- 14Y t 



m 2 

tL 


~ m\ yl 




~ m\ ^1 


A t 


~ ^1 


m t 





yk)~ 2mi /2 



(1.25) 



1 1 Several unification scales may arise as the gauge coupling unification scale and the string scale only few orders below the 
Planck scale. Numerically the variation of the scale Mx between these physical scales leads to a small effect only. 
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Figure 1.4: Some relevant masses in the approximate mSUGRA scenario for Aq = 300 GeV, tan /3 = 
4, fj, > 0; rriQ and m 1( / 2 are varied between 50 and 400 GeV. 



The IR fixed point of the top mass is Y" 4 IR ~ 8«3/9 and Y t /Y^ R varies from 0.75 to 1 dependent on 
tan (3, becoming unity for tan (3 = 1. In this limit the universal scalar mass does not influence the lighter 
right handed stop mass. If the doublet soft breaking mass is larger than the right handed soft breaking 
mass, the t\, defined as the light stop, will be mostly right-handed and the angle will prefer values around 
tt/2. 

The higgsino mass parameter in this limit will be given as 

r) 777 f O 1 O ^£ 

H + = -m - 2 m i/2 + terms including (1-26) 

The analyses in the following chapters are carried out using this approximate mSUGRA renormalization 
group solution^. If not explicitly stated otherwise we will vary the high-scale parameters around one 
central point: 



%/2 


=150 GeV 


m 


=100 GeV 


A 


=300 GeV 


tan (3 


=4 


n >0 


/'- 


=277 GeV 


M 2 


=122 GeV 


A t 


=355 GeV 










=55 GeV 




=103 GeV 




=283 GeV 




=309 GeV 




m xt 


=100 GeV 
=401 GeV 


m xt 
rriq 


=307 GeV 
=352 GeV 


ti 


=198 GeV 




=427 GeV 


sin(2^) = - 



0.97 (1.27) 

In Fig. |L4j some relevant low energy mass parameters are given as a function of mo and m^j to illustrate 
the qualitative behavior described above. Typical features are the large mass difference between the stop 
mass eigenstates, nearly independent of the value Aq, and the clustered neutralino masses, where the 



12 This is implemented in the initialization routine of SPYHTIA [B4]. Some comments concerning the 5.7 version can be 
found in the bibliography. 
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two light states are gaugino-type and the two heavy states are higgsino-type. The latter results from 
the large value for \i in the mSUGRA scenario. The lightest Higgs mass in this scenario in the given 
approximation is larger than 100 GeV and will not be excluded by LEP2. 

1.4. Mass Spectrum and Experimental Limits 
Neutralinos and Charginos 

Searches for neutralinos and charginos have been carried out at the Tevatron ^ as well as at LEP [Q]. 
Due to low energy R parity conservation they can only be produced in pairs XiXp XiXj-> m ^ X.fx'j- 
If the lightest neutralino is the LSP, then the heavier particles have to decay via a cascade into the LSP 
However the two and three partem decay channels are strongly dependent on the mass spectrum: 

x] — > z*xl h*xl — i + rxl QQX°i 

X] — d, vv 

Xl — W^xl H+*xl ■ ■ ■ 

x]-^tv,vl ■■■ (1.28) 

The decay x°j ~^ QQ Wli l be dominant if kinematically allowed, but in a SUGRA inspired mass scenario 
this will be only the case for the two heavy neutralinos. Besides, the chargino can enter the neutralino 
decay chain via Xj XiH~,xfW~. One very promising final state for the mixed neutralino/chargino 
production is the trilepton event 

PP/PP ^Xtx° 2 ^ tvtx Uxi — M + $r (1 -29) 

where three charged leptons are present in the final state and the missing transverse energy $t is based on 



three invisible particles. The exclusion plot is given in Fig. L5 , The cross section for chargino/neutralino 
production times the branching ratio into the trilepton channel is given for different squark masses, the 
gluino mass is fixed by the neutralino/chargino mass and the gaugino mass unification. The mass limits 
for xj can be read off the axis, they vary between 60 and 80 GeV. 

S quarks and Gluinos 

The gluino will in general be assumed heavy, as suggested by SUSY-GUT scenarios. The experimental 



exclusion limits from the direct search for squarks and gluinos are given in the mass plane in Fig. 1.5. 
The absolute lower limit on the gluino mass is m- g > 180 GeV rf7h. The decay channels considered for 
the light-flavor squarks and for the gluinos are 

Q — ► QX°j, q'xj — >jets + ^ T H 

g — > qqx°j, q'qXj — >jets + ^ T H — 

9^ hi ~^b X p--- (1.30) 

The final state neutralino/chargino decays via a cascade to the lightest neutralino, which is assumed to 
be the LSP. Products in this decay chain are denoted by the dots. If it is not kinematically forbidden, the 
gluino can first decay into a squark and a quark, and vice versa. This leads to one more jet in the final 
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Figure 1.5: Left: The CDF limits on a • BR/or the X2X1 production from the search for trilepton events, 



eq.( 1.29 ); Right: The CDF limits on the squark and gluino mass including the (jets+^x) and the like- 
sign lepton signal for the gluino pair production /|^/. In part of the parameter space the NLO cross 
sections have been used /[|/. 

state. The stop decay channel of the gluino leads to a higher multiplicity of Standard model particles and 
bottom jets. A typical signature for the Majorana gluinos arises from the decay via a chargino. Since the 
gluino is a singlet under the electro-weak gauge group, it decays to xf and xj with the same probability, 
leading to like-sign leptons in the final state of gluino pair production. A considerable Standard Model 
background is not present for this signature. 

Since supergravity inspired SUSY-GUT relations are used for the experimental search at hadron 



colliders, there are no strong limits on the squark mass if the gluino mass exceeds 550 GeV, see Fig. 1.5. 



The supergravity inspired GUT scenarios as described in section |L3| do not allow for a gluino mass being 
much larger than the light-flavor squark mass. In this region of the (rriq — nig) plane only the general 
unification of the gaugino masses can be kept. The mass of the lightest neutralino, assumed to be the 
LSP, grows with the gluino mass and becomes large enough for the squark to decay into an LSP almost 
at rest. The missing transverse momentum would then become too small to be measured. 

The limits on the neutralino/chargino mass from the search at LEP could be translated into limits 



on the gluino mass, using the gauge coupling unification eq.(1.22). Those are much stricter than the 
Tevatron limits but model dependent. 

Stops 

The limits on the stop mass arise from a search for stop pairs decaying into t\ — > c\\ and are therefore 
strongly dependent on the mass of the lightest neutralino. For a light stop mass this decay mode will be 
dominant. In this mass regime the light stop can be produced at LEP e + e~ — > t\t\, Fig. [L6| [gj. The 
production cross section depends on the mixing angle, arising from the t\i\Z coupling, and thereby also 
the mass bound. As will be shown in chapter ||, the hadroproduction cross section is independent of the 



18 



CHAPTER 1. SUPERS YMMETRY 



95% Exclusion Limit on t f — >c%° 




t, (GeV) 



Figure 1.6: Mass limits from the t\ pair production at DO and LEP followed by the decay t\ — ► cxi- The 
dependence on the mixing angle enters through the coupling t\t\Z at LEP. 



mixing angle, and both analyses, at LEP and at the Tevatron, yield a mass bound on the lightest stop 
around 79 GeV []9[ ^]. However, these limits are only valid as long as the \\ i s light enough and the 
decay channel t\ — ► c\\ is dominant. Additional limits arising from the search for the decay t — > iiXi 
are strongly dependent on the branching ratio of this decay mode and therefore weaker than those from 
the direct search. The different stop decay modes are described in chapter ||, and the direct search at 
hadron colliders is investigated in chapter [|. 



1.5. Regularization and Supersymmetric Ward Identities 
Dimensional Regularization and Reduction 

The MS renormalization scheme is by definition related to the regularization of infrared and ultravio- 
let divergences in dimensional regularization (DREG) [|35|]. This regularization scheme respects gauge 
symmetry and therefore the gauge symmetry Ward identities^]. It is less well-suited for supersymmet- 
ric theories, since all Lorentz indices are evaluated in n dimensions, whereas the spinors are still four 
dimensional. This leads to a mismatch between the degrees of freedom carried e.g. by a physical gluon 
(n — 2) and a gluino (2). A modified dimensional reduction scheme (DRED) has been introduced to cope 



with this problem [|37|]. The number of space-time dimensions is compactified from four to n dimen- 
sions, leaving the number of gauge fields invariant i.e. the gauge fields carry the n dimensional Lorentz 
indices. The remaining (4 — n) dimensions form the e scalars. These particles render the 7 algebra four 
dimensional. The gauge bosons and the gauginos carry the same number of 4 degrees of freedom. The 
DRED scheme will be used to illustrate the modified MS scheme. Except for the unsolved problem of 



mass factorization in DRED [38, m it can be shown that both dimension based schemes are consistent 



for calculations in the framework of supersymmetric gauge theories. 



13 We will not focus on the problem of the chiral projector matrix 75, since consistent schemes have been developed [|3^, |3^] to 
deal with 7 traces in n dimensions. In one-loop order a naive scheme can be used, however for neutralino/chargino production 
it has explicitly been checked that the ambiguous scheme dependent terms do not contribute [p5|] . 
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Starting with a Lagrangean £[W^,\ a ,D] for a non-abelian supersymmetric gauge theory in the 
Wess-Zumino gauge one can show that the supersymmetric variation 5s of the Lagrangean only vanishes 
in the limit of (n — ► 4) dimensions, up to a total derivative [ ]3^ ] 

5 s C[W«,\ a ,D] (1.31) 

The component fields W indicate the gauge fields, and A the [Majorana] gauginos. This leads to the 
Ward identity including the ghost and gauge fixing term Lq, where in n dimensions the variation 5s£ 
has to be kept. 

= ( J d n x [J" 6 S W„ + ]5 S X + 3d S S D + 5 S C G + 5 S C] ) 

(X}= J d{W^}d{X}d{D}X e iId n xlC+C G +J^+]\+j D D] (L32) 



Although DREG and thereby the MS scheme cannot be shown being inconsistent with supersym- 
metry, they do not respect supersymmetry on the level of naively used Feynman rules. The problem is 
similar to applying DRED to gauge theories: Evanescent couplings renormalize in a manner different 
from the physical couplings . In NLO-DREG this results in a finite renormalization of Feynman diagrams 
which restores supersymmetry explicitly. At higher orders these additional counter terms even include 
poles in e. 



Finite Renormalization 



Explicit calculations show that Green's functions calculated from the MSSM Lagrangean using di- 
mensional regularization may not respect supersymmetry. The supersymmetry transformation mirrors 
e.g. the gauge coupling g{qqg) to the gauge coupling g(qqg) and the Yukawa coupling g(qqg). In reg- 
ularization schemes which respect supersymmetry, like dimensional reduction^, the supersymmetric 
limits of these couplings are identical in any order perturbation theory. In DREG the supersymmetric 



limit of the Yukawa coupling differs from the gauge couplings at one loop level [ 4C ] 



9 = 9 



1 + 



32vr 2 



\C{G) 



(1.33) 



The Casimir invariants C are defined for the Dirac fermions in the fundamental (r), and for the gauge 
boson and the Majorana gauge fermions in the adjoint (G) representation]^]. This difference has to be 
compensated to render the calculation supersymmetric. Since the Standard Model quark-gluon coupling 
g(qqg) is by definition the measured quantity, the Yukawa coupling will be shifted g — > g in the ex- 
pression for the final observable. This finite shift is not a finite field theoretical renormalization of any 
measured parameter and it is not only present for gauge vs. Yukawa couplings. It is an artifact arising 
from the supersymmetry violation of naive dimensional regularization. 



14 The difference between DREG and DRED are e terms arising from DREG Dirac traces including gauge fields. They 
combine with a pol e 1/e in a scalar integral, leading to a finite contribution. These terms are exactly those leading to the 
difference e.g. in eq.(1.34). 

15 The SU(3) coupling qqg yields C(r) = C F and C{G) = C a = N. 
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Supersymmetry relates the weak Higgs Yukawa coupling Y(qqh) to the vertices Y(qqh) and Y(gqh). 
The three couplings Y in the supersymmetric limit and calculated in DREG are not identical in NLO 



Y(qqh) = Y(qqh) 



1 + 



167T- 



:C(r) 



Y(qqh) 



1 + 



*2L 

32vr : 



■C(r) 



(1.34) 



where in the case of weak coupling only C(r) occurs. These two finite differences oc a s in couplings 
mediated by Gf have to be compensated to make dimensional regularization compatible with supersym- 
metry. 

The usual parameterization of the Yukawa coupling constant is Y = mg, where g is defined in the 
MS scheme and m is the pole mass, i.e. renormalized in the on-shell scheme. However, the pole mass 



has to be calculated in the DREG scheme, and the mass appearing in the different couplings eq . d 1 . 34j ) is 
— in the supersymmetric limit — only numerically the same. In fact, the scalar mass set to m and the 
fermion mass set to m behave differently in next-to-leading order, since the counter term for the scalar 
and the fermion on-shell mass in DREG is not the same. 



m„ 



1 + 



9 



16vr 2 



C(r) 



ma 



(1.35) 



This behavior breaks supersymmetry explicitly and has therefore be removed. The mass shift is respon- 
sible for the difference between Y(qqh) and Y(qqh), and it renders the difference to Y{qqh) compatible 
with the general difference between the gauge and Yukawa coupling as given in eq.( |l.33| ): 



g(qqh) = g{qqh) = g(qqh) 



1 + 



32vr 2 



C(r) 



(1.36) 



The observable coupling is again defined as in the Standard Model value Y(qqh). 



2. Production of Neutralinos and Charginos 



2.1. Born Cross Sections 
Partonic Cross Sections 

Neutralinos and Charginos can be produced at hadron colliders in several combinations, all starting from 
a pure quark incoming state 

QQ — * XiX j 
QQ — ► xfxj 
ud — * xtx°j 

du — XiXj (2-1) 
The first two processes are possible for a general quark-antiquark pair. For the latter, charge conservation 



requires u and d type quarks in the initial state [41] 



Two generic Born Feynman diagrams contribute [Fig. [Q| ]: an s channel gauge boson (7, W,Z) 
Drell-Yan like and (t, u) channel squark exchange diagrams. Two final state neutralinos are produced 
by the first diagram purely as higgsino-type. Final state charginos can couple to the s channel gauge 
boson as gauginos and as higgsinos. For mixed neutralino/chargino production the s channel diagram 
contributes to all current eigenstates as well. In the given approximation of a trivial squark CKM matrix, 
the t, it channel squark couples flavor conserving to the incoming quark and will therefore be regarded 
as light-flavored; the incoming quark originates in the parton density of the proton, and will consistently 
be assumed massless. This makes the higgsino Yukawa coupling vanish for all possible final states. For 
the gaugino-like charginos this coupling also vanishes in case of q R , since the qq\ coupling respects the 
helicity eigenstates. 

The LO partonic cross section a, which is proportional to the matrix element squared in the limit of 
(n — > 4) dimensions can for all possible final states be written as [The n dimensional Born cross section 




f q 

-e WW 



Figure 2.1: Generic Born diagrams for neutralino/chargino production 
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is required for the NLO contribution.] 



da. 



dt 



2ira 2 
Ns 2 



2mm i (A Ll A* Lj C T2 + A Ri A Rj C T4 ) + ^ {A Ll A* Lj C T1 + A Rl A* Rj C T3 ) 



tn 



(A c Lj A c * Ll C T1 + A c Rj A c * Ri C T3 ) - ^ (A c Lj A c * Li C T2 + A^A^CW) 



+ 



t{tj + //,•//,• "';"' / „ 
5 ^51 H ^52 + 



2{u x u 2 - tit 2 ) 
s(s-M 2 ) 



Cs3 



8smirrij 

tqUg 
8titj 



(A Lj Al t A c Lj A c * Lt + A Rj A Ri A c Rj A c * R i) 



+ -2T {ALiAliALjAlj + AjaA^AjyA^-) 
+ (A^Cfc + Ajy4fcC$ 4 ) + ^ (A Lj A^ + A Rj A* m C} 3 ) 



8smimj 



tnU, 



qu-q 



{A Li A%A c Li A c % + A ra A^A c ffi A%.) 



+ 



8td 



i '.i 



+ 



2mim 



[A c Ll A c * Ll A c L3 A c * Lj + A c Ri A c * Ri A c R] A%) 



2.1. I ~£> n 



Ur, 



(A c Ll A c * L3 C* T1 + A c Ri A%C* T3 ) + -U- (A c Ll A c * Lj C* T2 + A c m A c * Rj C* TA ) 

SUfj 



where 



tg = t- mq 



ur, = u — nip, 



(2.2) 



M is the W or Z mass of the s channel gauge boson. The coupling parameters A correspond to the 
t, u channel couplings for the outgoing particles i, j and are defined in Tab. |A.5| . The charge conjugate 
coupling A c is identical to A for the neutralinos. In the chargino case Aj is the coupling for outgoing 
Xj containing the mixing matrix U, and A c j for an outgoing \j containing V. The typical couplings 



Cs, Ct follow from the Feynman rules Fig. |27T 

X rS 



Csi = X c 



2(s — M 2 ) 



Re [{£ + r)(L + R)\ + 



4(s-M 2 ) 2 

2 



Cs2 = X£. 



X '" Re[(e + r)(L + R)} + 



2(s — M 2 ) 



Cs3 
Cti 



— Re \(e-r)(L-R)] + , - ^ 
4 LV A JS 8(s-M 2 ) 



X c + 



Ct3 = X c + 

x r = 



s-M 2 
s 



IR 
rL 



Ct2 = X c + 



2(s-M 2 ) 2 

2_| r |2)(| L |2_| i2 |2 
S 



(|^| 2 + |r| 2 )(|L| 2 + |i?| 2 ) 



(|£| 2 + |r| 2 ) Re [|L||iJ|] 



s-M 2 
Q only for xJXj 



Ct4 = X c + 



s-M 2 
s 



M 2 



rR 



(2.3) 



The gauge boson-quark couplings r, £ are given in Tab. 



the neutralino-chargino couplings in 



Tab. A.4. Final state charginos require one subtlety in the matrix elements: either the t or the u chan- 
nel diagrams contribute to the amplitude with a fixed quark flavor, except for the pure neutralino case. 
For two final state charginos, A only couples to u type, A c to d type quarks. In the mixed production 
processes the couplings A and A c have to be arranged making use of charge conservation. 
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The factor of C53 in the Born cross section eq.([2T2j) originates from the contraction of two CP odd 
Dirac traces Tr(757^i7^2)Tr(757 / ^i7 l/ ^2), where the definition of the momenta is given in ap- 



pendix B.l. Using a naive 75 scheme, this term cannot be fixed consistently. We therefore keep this 
kind of structure in the Born, real gluon and virtual gluon contributions. The different choices for the 75 
scheme result in 0(e) corrections and do not contribute to the final expression, since the corresponding 



diagrams are finite. The calculation performed in the consistent 't Hooft-Veltman scheme [35] agrees 
with the naive calculation. 

Hadronic Cross Section 

The hadronic cross section for pp/pp collisions is given by a convolution of the partonic cross section 
with the parton densities for the quarks in the proton, e.g. for two hadrons H\H 2 



a(S,Q 2 ) = [ dx i [dx j [f^(x l ,Q 2 )f^(x 1 ,Q 2 ) + (H 1 ^H 2 )\a i:j (x l x J S,Q 2 ) 

partons ij T ° x i 

(2.4) 



where k\ and k 2 are the incoming parton momenta, S = (k\ + k 2 ) 2 is the hadronic cm energy; rrij are 

the masses of the final state particles, and to = (mi + m 2 ) 2 / 'S is the kinematical limit, ff 3 are the 
parton densities, forming the convoluted hadronic luminosity 



a(S,Q 2 )= £ | 1 d r ^( T ,Q2 ) ^.. (r5) Q2^ 

partons ij T ° 



^(r, Q 2 ) = [fi ® f d ] (r, Q 2 ) + [fj M (r, Q 2 ) 

[f ® 9] (r, Q 2 ) = f 1 — f(x, Q 2 ) 9 (-,Q 2 ) (2.5) 

where the hadrons Hi , H 2 are implicitly fixed by the order of the convolution of the parton densities. For 
identical incoming gluons a factor 1/2 has to be incorporated. 

2.2. Next-to-leading Order Cross Sections 
2.2.1. Virtual and Real Gluon Emission 

The NLO cross section includes the radiation of real quarks and gluons and virtual gluons and gluinos. 



The generic diagrams are given in Fig. 2.2 for the qq incoming state. The additional qg and gq diagrams 



are obtained by crossing one quark to the final and the gluon to the initial state. The virtual contributions 



are regularized by dimensional regularization. Therefore a finite shift of the couplings eq.( |1.34| ) has to 
be applied to restore supersymmetry. The divergences appeal - as poles in e, as shown in appendix B.3. 
The UV poles require renormalization; the only parameter in the Born term eq.( |2.2| ) which undergo 
the renormalization procedure is the squark mass, defined as the pole mass, i.e. in the on-shell scheme. 
The soft gluon poles cancel with the real gluon emission. The phase space integration for the real 
gluon emission is given in appendix |B.1| . These matrix elements have been computed using phase space 
subtraction, i.e. the additional gluon phase space is integrated numerically. After subtracting the dipole 
terms the remaining divergences are of collinear type and removed by mass factorization, appearing in 
the subtraction term, see appendix 
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Figure 2.2: Generic NLO diagrams for neutralino/chargino production, the self energy contributions are 
not shown. 



2.2.2. Mass Factorization 

The partem densities eq.([2~3|) form observable structure functions [e.g. F 2 ], which contain divergences in 
next-to-leading order QCD p4|]. These divergences arise from the collinear radiation of gluons and have 
a universal structure which is fixed by the Q 2 evolution. They have to be absorbed into the definition of 
the parton densities to render the physical structure function finite. In analogy to a UV renormalization 
procedure it is possible to absorb additional finite parts into the re-definition. The minimal set is the MS 
scheme, and it leaves the next-to-leading order contribution to the measured structure function with a 
non-zero finite term. This minimal choice respects the required sum rules naively. 

Due to the factorization theorem, the universal form of the partonic cross section in the collinear limit 
is independent of the order of perturbation theory. 



2 d?°~ij 

3 — 

dt2dsi 



l dx> f'd^ 2 2 U^jjt 

Xi Jo Xj V dt 2 ds 4 



Sij 6(1 



!_ f^9l 

2ttc T(l - 2e) V Q 2 F 



r(i 



Pij(x) 



(2.6) 



Tij is called splitting function and describes the splitting of a parton i to a parton j in the collinear limit. 
It is evaluated perturbatively and consists of the trivial LO term and a divergent NLO contribution. The 
appearance of the Altarelli-Parisi kernels Pij fixes the Q 2 evolution, they are given in eq.(B.23). Other 



non-minimal schemes lead to a finite renormalization — > + fy. The reduced cross section <r is 
finite and, as well as the splitting function, depends on the factorization scale Qp. This scale dependence 
should flatten after adding higher order perturbative contributions, since it is a perturbative artifact. 
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Figure 2.3: Feynman diagrams for crossed channel production of neutralinos/charginos including on- 
shell intermediate states, which have to be subtracted. 



The renormalization of the parton densities has to cancel the remaining collinear poles in the matrix 
elements and leave the final expression finite. The counter term which has to be added to the bare cross 
section to obtain the reduced one in the MS scheme can be read off eq.(j2li|) 



2 d 2 ^ F = a s r(l-e) /4vrQ 2 y f 1 dx 
S dt 2 ds 4 2vrer(l-2e) \ Ql J J x 



2.2.3. On-Shell Subtraction 



dt 2 ds4 j V dt 2 ds4 

xkr 



i/m 



(2.7) 



Apart from the UV and IR divergences another kind of divergences can occur, due to on-shell inter- 
mediate particles. After crossing the NLO production matrix elements, different incoming states may 
contribute to the (xx+jet) inclusive final state 

qg — * XiXjQ 

gq — * XiXjq (2-8) 

As depicted in Fig. [2.3[ , these can proceed via an on-shell squark. A natural way of solving the prob- 
lem would be introducing finite widths for all particles under consideration. However, a finite squark 
width would spoil gauge invariance. In addition, it would yield a strong dependence of the next-to- 
leading order production cross section on the physical widths of intermediate states. This dependence 
would only vanish after including the decays into the calculation. Therefore we instead differentiate 
between off-shell and on-shell particle contributions, the latter regarded as final states in the set of super- 
symmetric production cross sections. 

Considering an analysis of all production processes for two MSSM particles at hadron colliders this 
differentiation removes a double counting of the on-shell contributions of the squark, as it would occur 
in the case of general finite widths: 

gq — > q*Xi — ► qXjXi neutralino pair production 

gq —> q~Xi ' BR(g — ► q\j) squark neutralino production (2.9) 
The on-shell squark contribution is subtracted from the crossed XiXj production matrix element, leaving 



it as a contribution to direct \Q. production, eq.(|2.9|). The off-shell contribution is kept for the first 



of the processes under consideration. To distinguish these contributions numerically, one regularizes the 
possibly divergent propagator by introducing the Breit-Wigner propagator (p 2 —m 2 ) — ► (p 2 — m 2 +imT). 
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Since this width can be regarded not as a physical property of the final state particle, but as a mathematical 
cut-off, the matrix element can be evaluated in the narrow width approximation, regarding the final state 
particles as quasi-stable. 

Assuming an on-shell divergence in the variable M 2 , the hard production cross section in the narrow 
width approximation reads 

da , _ . mgTg/ir . 

^ = a( 9g ^ qXl ) (M2 _ ^ )2 - BR (q -+ qXj ) + O 

— > * (59 - BR (g -> ,5(M 2 - m|) + ^^^2 J ( 2 - 10 ) 

In case of the neutralino/chargino production M 2 = S3 + m\ and M 2 = 54 + m 2 , are relevant for 
the on-shell squarks, the extended set of Mandelstam variables is defined in appendix |B|. The leading 
divergence is subtracted from the crossed channel matrix element, as described before. The complete 
crossed channel matrix element can be written as \M\ 2 = /(M 2 )/[(M 2 — m~) 2 + m|r|]; then the 
subtraction for an intermediate squark is defined as 

f(M 2 ) f{m 2 ~) 

~7~rrK 0-7; 9—n — ttt^ oTq tt^j @(s — (rria + rrii) )Q(m d — mj) (2.11) 

(M 2 - m 2 ) 2 + m?r? (M 2 - m 2 ) 2 + m\Y\ q 1 V q 3 

Since an over-all factor S(M 2 — m~) is absent in the subtracted term, the Breit-Wigner propagator has 
to be integrated over the phase space variable M 2 . The matrix element, including the remaining phase 
space integration is evaluated for M 2 = m 2 -. 
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Figure 2.4: The renormalization/factorization scale dependence of the total cross section for x\X\ pro- 
duction at the upgraded Tevatron and the LHC. There is no maximum in either of the next-to-leading 
order curves, and the LO and NLO do not meet for a scale around the average mass. The SUSY scenario 



determining the masses is given in eq.( 1.27). 
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The remaining non-leading divergences, arising from interference between finite and divergent Feyn- 
man diagrams, are integrable and well-defined using a principal-value integration. Numerically this 
principal value can be implemented by introducing a small imaginary part (M 2 — > M 2 — ie). Since the 
matrix element squared may contain subtractions in more than one variable this imaginary part may lead 
to finite contributions and has therefore to be taken into account. 



2.3. Results 



Scale Dependence 



Since the leading order hadro-production cross section for neutralinos and charginos does not contain 
the QCD coupling constant, it only depends on the factorization scale through the parton densities. This 
renders the leading order scale dependence smaller than ~ 30%. The variation of the cross section with 
the scale is therefore not a good measure for the theoretical uncertainty. In next-to-leading order, this 
factorization scale dependence becomes weaker; however, an additional dependence on the renormaliza- 
tion scale arises. For p, = pLp = p,R this yields a generally weak scale dependence of < 20% at the 
upgraded Tevatron and < 5% at the LHC. As can be seen from the leading order curves in Fig. 2.4, the 
combination of factorization and renormalization scale dependence leads to a different behavior at the 
Tevatron and at the LHC, due to different momentum fractions x contributing; in contrast to the strong 



coupling induced processes a maximum cross section for some small scale does not occur, Fig. 2.4 
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Table 2.1: A complete set of K factors for neutralino and chargino production at the LHC. The masses 
are chosen according to the default SUGRA inspired scenario, eq.{ 1.2^ . The renormalization and fac- 
torization scales are set to the average final state mass. Although the K factors are of a similar size 
1.3 •• • 1.5 for each diagram contributing, large cancelations lead to huge corrections for the scenario 
under consideration. 



28 



CHAPTER 2. PRODUCTION OF NEUTRALINOS AND CHARGINOS 



Numerical Results 

The production of neutralinos and charginos can be probed at the upgraded Tevatron, a pp collider with a 
center-of-mass energy of 2 TeV, and at the future LHC, app collider with an energy of 14 TeV. The cross 
section for several combinations of light neutralinos and charginos, which turn out to be gaugino-like in 



the considered scenario, are given in Fig. |2.5[ The size of the cross sections strongly depends on the 
mixing matrix elements associated with the different couplings. This yields e.g. a larger cross section for 
X2X1 pairs compared to XiXi production. In general, the processes containing no final state chargino 
are suppressed, independent of the masses, which are almost the same for an(1 X\- Whereas the cross 
section for the production of positively and negatively charged mixed pairs are identical at the Tevatron, 
they differ significantly at the LHC, due to non-symmetric parton luminosities. The dependence on 
SUSY masses and parameters, which are not contained in the leading order cross section, like the gluino 
mass, is weak in next-to-leading order. The virtual corrections are generically small [< 10%] compared 
to the real gluon emission; however, they are not universal and even do not have a unique sign for the 
different gaugino and higgsino-type outgoing particles. 

The next-to-leading order K factor is consistently defined as K = <jnlo/&lo- It is dominated by the 
gluon emission off the incoming partons and therefore similar for all considered processes and a constant 



function of the masses, Tab. [2.1[ Although the real gluon corrections to any diagram contributing to the 
production process are of the order 1.3 • ■ • 1.5, large cancelations give rise to huge K factors. The same 
effect occurs for the virtual corrections, which grow up to 50% e.g. for the X2X2 or the X4X1 channel. 
Varying the common gaugino mass m^j reduces the K factor to values expected by regarding the other 
channels. 

With an integrated luminosity of J C = 20fb _1 in run II, the upgraded Tevatron will have a maximal 
reach for the mass of the produced particles when probing the X2X1 channel. For masses smaller than 
150 GeV, 10 3 to 10 5 events could there be accumulated. Although the XiXi cross section is compatible 
with the mixed neutralino/chargino channel for a fixed value of the common gaugino mass, the particle 
masses, which can be probed, stay below 80 GeV in the considered SUGRA inspired scenario. The 
same holds for the LHC, where for typical masses of the X\ an(1 X2 below 300 GeV and an integrated 
luminosity of f C = 300fb _1 a sample of 10 4 to 10 6 events can be accumulated. In the given scenario 
the higgsino type neutralinos and charginos are strongly suppressed compared with the lighter gauginos. 
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Figure 2.5: Some total cross sections for pairs of neutralinos and charginos at the Tevatron and at the 
LHC as a function of mi/ 2 - The different masses of the particles involved are given on secondary axis. 
The strongly suppressed heavy higgsino cross sections are not given. The and x\x{ cross sections 
are different due to the non-symmetric parton luminosities. 



3. Scalar Top Quark Decays 



Scalar top quarks can decay into two or three on-shell particles via the strong or electroweak cou- 
pling [Qj. The possible two body decays are — kinematically allowed for an increasing stop mass 
in typical mass scenarios: 



W + b x 



bx+, txl [Zt u hh, W + b, H + b, tg] 



i- 



H + bx\ 



(3-1) 



The channels in brackets are possible only for the heavier stop, since the t\ is assumed to be the lightest 
scalar quark. The decay into a charm jet is induced by a one-loop amplitude, and will therefore be 
suppressed, if any other tree-level two or three body decay channel is open. In the intermediate mass 
range, when the b\f channel is still closed, the three particle decay into Wbxi is dominant. For a heavy 
I2 the strong decay mode including a final state gluino will be the leading one, as will be shown later in 
this chapter. 



3.1. Strong Decays 
3.1.1. Born Decay Widths 

Since the Yukawa qqg couplings are flavor diagonal, any decay involving a scalar top quark 



v 3 

9 



t + g 

t + tj and c.c. 



m^. > nit + m g 
nig > nit + nif. 



(3.2) 



includes a top quark in the final state, i.e. the strong decays will only be possible for large mass scenarios. 



For the light stop t\ the weak decays in eq.( |3.1[ ) will be the only kinematically allowed. 

The calculation including the stop mixing and a massive top quark is a generalization of the light- 
flavor decay width [fl||]. To lowest order the partial widths for the stop and gluino decay, eq.(32), are 
given by by[] 



r( ^ 2 t§) = iv^ Al/2 K, 2 '™?'^) 
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(3.3) 



The different factors in front are due to the color and spin averaging of the decaying particle, and the 
crossing of a fermion line. Interchanging t\ and ti in the two leading-order decay widths corresponds to 



the symmetry operation Vyi in the Lagrangean, as described in section 1.2.4. 



'A(a;, y, z) = x 2 + y 2 + z 2 — 2(xy + xz + yz) 
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Figure 3.1: (a) Born diagrams for stop and gluino decays; (b) vertex corrections; (c) real gluon emission. 
The correction to the gluino decay can be obtained by crossing the diagrams 



3.1.2. Next-to-leading Order SUSY-QCD Corrections 
Massive Gluon Emission 

The NLO corrections [|ll|] include the emission of an on-shell gluon, Fig. |3.lfc . This gluon leads to IR 
singularities which are regularized using a small gluon mass A, subsequently appearing in logarithms 
log A 2 . The massive gluon scheme breaks gauge invariance for the non-abelian SU(3) symmetry. Hence 
the scheme has to be extended by new counter terms if a non-abelian contribution arises from a three or 
four gluon vertex, otherwise the SU(3) Ward identities would not be satisfied anymore. This is not the 



case for the stop decays Fig. p.l[ The gluon behaves like a photon and its mass can be regarded as a 
mathematical cut-off parameter. After integration over the whole phase space the small mass parameter 
drops out and yields a finite sum of virtual and real gluon matrix elements. However, these massive 
gluon matrix elements must not be interpreted as exclusive cross sections, since gauge invariance is only 
restored for inclusive observables, i.e. the gluon integrated out. 

In the considered process the logarithms of the gluon mass arise from the integration over the soft 



and collinear divergent three particle phase space, eq.(B.7). The same kind of logarithms enter through 



the virtual gluon contributions, e.g. the scalar three point function eq.( |B.44] ) and cancel analytically. 



Virtual Corrections 

The virtual gluon corrections, including self energy diagrams for all external particles and vertex cor- 
rections Fig. [3.1| s, are also regularized using the massive gluon scheme. The additional UV divergences 
have to be regularized dimensionally. The poles 1/e are absorbed into the renormalization of the masses, 
the strong coupling, and the mixing angle, which are the parameters appearing in the Born decay width 



eq.(J33j). The counter terms for mass renormalization in the on-shell scheme and the renormalization of 



a s in MS can be found in appendix B.4. The mixing angle is renormalized by introducing the running 



mixing angle and absorbing the mixing stop self energy contributions. This scheme restores the {t\ «-> £2) 
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H/m(g sin(2e) 



Figure 3.2: Left: Renormalization scale dependence of the gluino decay width for the same SUGRA 
inspired scenario as chosen in Fig. [Q| ; mi = 449GeV,mj 2 = 847 GeV, sin(2#) = —0.59, and 
m g = 637 GeV. The renormalization scale is varied as a fraction of the mass of the decaying gluino; 
Right: mixing angle dependence of this gluino decay, where the stop mixing has been varied over the 
whole range, independent of the other low energy parameters. 



symmetry Vn in NLO. The dependence on the mixing angle in NLO can be described by a constant K 
factor, Fig. |3.2| , possibly large contributions from the gluino-top loop are absorbed into the definition of 
the mixing angle. Renormalizing the strong coupling in the MS scheme breaks supersymmetry; adding 



a finite counter term, derived in eq.(1.33), restores supersymmetry. 



The Born decay widths are proportional to A 1 / 2 , i.e. the relative momentum of the produced particles. 
One of the vertex correction diagrams is constructed by exchanging a virtual gluon between outgoing 
color charged particles. Near threshold the exchange of a gluon between two slowly moving particles 
picks up a factor A -1 / 2 , the Coulomb singularity, which cancels against the phase space suppression 
factor in the virtual correction matrix element. The NLO decay width therefore does not vanish at 
threshold. The narrow divergence can be removed by resummation of the contributions near threshold. 
Moreover, the screening due to a non-zero life time of the final state particles reduces the Coulomb effect 
considerably. 

The complete analytical expression for the stop decay width is given in appendix |C]. The numerical 
results are shown together with the weak decays in Fig. 



3.2. Weak Decays 

The possible weak decay modes including a stop will be dominant once the strong channels are kinemat- 
ically forbidden. Although this region is not preferred by the mSUGRA scenario even the crossed top 
decay could be possible, which leads to experimental limits on the branching ratio of this decay mode 
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and thereby on the masses involved |7|] 











t - 





rrif. > rrib + 
[mt > m h +m^] 



The Born decay width for the t\ decay to a neutralino reads^ 



2a 



A 1/2/ 2 2 2 \ 

■A / (m fi ,m t ,m^.) 



2 2 

m,- — m f 



(3.4) 



m~. j (Cf + C&) + Amtm^CLCR 



C L = A L cos 6> + B L sin 6* 
Cx = Br cos — Ar sin 



(3.5) 



The couplings A and i? are given in Tab. |A.5| for the neutralino involved. The t-i decay can be derived 
using the V\i operation. The decay channel producing a bottom quark and a chargino can be read off 
using Tab. |A.5| by setting the mixing angle to zero, as long as sbottom mixing is neglected. The NLO 
calculation is performed exactly as for the strong decay channel, whereby some virtual and real correction 



diagrams in Fig. 3.1 vanish for a Majorana particle without color charge. Again the finite shift eq.(1.34) 
has to be added to the weak coupling vertex, no matter if a gaugino or a higgsino is involved. 



3.3. Results 



In the calculation the renormalization scale of the process is fixed to the mass of the decaying particle. 
Since the scale dependence should vanish after adding all orders of perturbation theory one expects the 



variation of the width with the scale to be weaker in NLO than in LO. This is shown in Fig. for the 
strong coupling gluino decay. 



The numerical results for the strong decay channels can be seen in Fig. 3.3. Assuming for illustration 
a SUGRA inspired mass spectrum the light stop t\ can decay only via the weakly interacting channels. 
The strong decays are possible for the gluino and for the heavy stop. With increasing my 2 tne gluino 
becomes heavier compared with the stop masses, i.e. the decay into the gluino vanishes and the gluino 
decay channel g — ► t\t + t\t opens. A kink in the NLO ti decay widths occurs at the production 
threshold g — > tb\, where the gluino self energy exhibits a large discontinuity. It can be smoothed out by 
introducing a finite width of the gluino. The Coulomb singularity is present in both of the strong decay 
channels. However, it can be seen only in the stop decay, since the kink in the gluino self energy and the 
Coulombic vertex contribution to the NLO decay width cancel each other numerically near threshold. 
Since each of the large contributions is narrow, the phenomenological consequences are negligible. 

The large difference in the size of the virtual corrections between the stop and the gluino decay is due 



and arise through the analytical continuation 



to 7r 2 terms which are determined by the sign of (mj — nig 
of the matrix element squared into the different parameter regions, Fig. 3.3. For the gluino decay they 
give rise to destructive interference effects of the different color structures, and render the over-all NLO 
collections small. The size and the sign of the NLO correction to the gluino decay depends on the masses 
involved. The K factor for the stop decay is always large and positive K = 1.35 ••• 1.9 decreasing far 



"This decay width has also been calculated in NLO by other groups [|29|]; we have analyzed it for the sake of comparison 
and to illustrate the running mixing angle. The three calculations are in agreement. 
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Figure 3.3: Left: SUSY-QCD corrections to the strong decays t2 — > gt and g — > tt\ + tt\ as a function of 
the common gaugino mass. The masses of the particles involved are labeled in the additional axis. The 
input parameters have been chosen as thq = 900 GeV, Aq = 900 GeV, tan (3 = 2.5, p, > 0, in order to 
see the possible structures of the curves. The kink in the NLO stop decay width results from the gluino self 
energy and could be smoothed by inserting a finite gluino width; Right: SUSY-QCD corrections to weak 
decays of the light t\. The input parameters are the usual mo = 100 GeV, Aq = 300 GeV, tan (3 = 
4, /j, > 0. Only the decay into the two lightest neutralinos and into the light chargino is possible in the 
mi/2 range considered. Note that the mass of the \t an d tne X2 are a l most identical in the SUGRA 
inspired scenario. The dashed line denotes the LO, the solid one the NLO results in both figures. 



above threshold, the K factor for the gluino decay is in general modest and tends to be smaller than one, 
K = 0.8---l. 



The weak decay widths of the light stop t\ are shown in Fig. 3.3. They are genetically suppressed 
compared to the strong decay widths, due to the coupling constant. This yields about one order of 
magnitude between the different contributions. Moreover the typical weak coupling factor includes 
mixing matrices of the neutralinos and charginos, which may lead to a further suppression. Given that 
the masses of the four neutralinos cluster for the higgsino type and for the lighter gaugino type mass 
eigenstates, even the decay width into the heavier neutralinos/charginos can exceed the width to the 
lighter one [11]. Since the top quark is heavy, the bxf decay mode is typically the first tree level two 
particle decay kinematically allowed. The neutralino channels open only for higher stop masses, but 
will then be of a comparable size. The NLO corrections exceed 15% for special choices of masses and 
parameters only 
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3.4. Heavy Neutralino Decay to Stops 



Heavy neutralinos will be produced at a future e + e linear collider [[K]]. In most supergravity inspired 



scenarios they are higgsino-like, and will therefore not decay into light-flavor quark jets. However, the 

or a light stop. The 
]5|), by crossing the 



large top Yukawa coupling may open the decay channel — ► t%t + t\t [j = 3, 4] for a light stop. The 



analytical expression for this decay can be obtained from the stop decay width, eq. 
stop and the neutralino. 
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l + cl) 



+ Am t m^C L C R 



(3.6) 



The couplings are defined as for the stop decay. The numerical result is shown in Fig. 3.4. Similar 



to the strong decay width the neutralino decay to two strongly interacting particles exhibits a Coulomb 
singularity, due to the exchange of a slowly moving gluon between the decay products. Especially for the 
steeply rising decay of the second heaviest neutralino, this Coulomb singularity is very narrow, as can 



be seen in Fig. 3.4. Unlike the gluino decay in section 3.1 the next-to-leading order corrections can be 
of the order 10% and are positive over the whole mass range. The clustering of the masses of the heavy 
higgsino-like neutralinos is a typical behavior in supergravity inspired scenarios, since the off-diagonal 
entries of the neutralino mass matrix are small compared to \x. 
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Figure 3.4: SUSY-QCD corrections to the weak decays of heavy neutralinos to a stop x*j — ► tt\ +tt\ as a 
function of the common gaugino mass. The masses of the particles involved are labeled in the additional 
axis. The input parameters have been chosen as toq = 450 GeV, Aq = 600GeV,tan/3 = 4, fj, > 0, 
in order to see the Coulomb singularity at threshold. The dashed line denotes the LO, the solid one the 
NLO results. 



4. Production of Scalar Top Quarks 



4.1. Diagonal Stop Pair Production 



4.1.1. Born Cross Section 



Diagonal pairs of stop particles can be produced at lowest order QCD in quark-antiquark annihilation 
and gluon-gluon fusion: 



qq 
99 



and 
and 



(4.1) 



Mixed pairs ?i?2 and i 2 ii cannot be produced in lowest order since the git and ggti vertices are diagonal 
in the chiral as well as in the mass basis of the squarks involved. The relevant diagrams for the reactions 



(4.1) are shown in Fig. 4.6 a. The corresponding cross sections for these partonic subprocesses may be 
written as 

a LO [qq^tt] = — (3 3 
s 27 



olo [gg 



it] 




(4.2) 



The invariant energy of the subprocess is denoted by y/s, the velocity by (3 = 
sections coincide with the corresponding expressions for light-flavor squarks 



1 — 4m?/s. The cross 



<5lo [<?<?' -> qq] 



2 



&Lo[gg 



2nqCt 2 s -K 



9(/x 4 + m|s) 



+ 



+ 8, 



qq 



4 V 
9 9s 



L 
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L = log((l -0 + 2fi 2 /s)/(l + (5 + 2fi 2 /s)) (4.3) 



in the limit of large gluino masses and for 2n,q = 1. The main difference between these two cross 
sections results from the flavor diagonal qqg coupling, which makes the t channel gluino contribution for 
the squark production eq.(4.3) vanish in case of stops. This yields a /3 3 dependence of the qq — ► ti cross 



section. As described in section |1.2.3| an additional factor 2riq arises for the mass degenerate light-flavor 
squark production. 
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Internal gluon propagators in the LO and the NLO calculation are evaluated in the Feynman gauge. 
External gluons are restricted to their physical degrees of freedom. The sum over the physical polariza- 
tions in the axial gauge is 



n^k v + k»n v n 2 k^k v 



transverse DOF 



o 



(nk) 
k„P^ 



(nky 



(4.4) 



with an arbitrary four vector n. In the final result this vector n drops out according to gauge invariance. 

-gt") and 



Using the Slavnov-Taylor identity ( |B.47D this is equivalent to using the polarization sum 
removing the momenta k± and k\ from the tensor matrix element M.^ v {k\, A^)- 



4.1.2. Next-to-leading Order Cross Section 

The incoming gluons in the virtual and real correction matrix elements are treated the same way as 
in the Born matrix elements. The Feynman diagrams for the virtual gluon correction are shown in 



Fig. 4.6 b,c. The masses are renormalized in the on-shell scheme, the coupling constant a s in the MS 
scheme. The renormalization is performed in such a way, that the heavy particles (top quarks, gluinos, 
squarks) decouple smoothly for scales smaller than their masses, as described in eq.( |B.49 ). Note that 
no vertex requiring a finite renormalization according to section L5 occurs in the Born term, again in 
contrast to the light-flavor squark production. 

The calculation of the gluon bremsstrahlung matrix element has been performed in the cut-off 



scheme, appendix 25 The Feynman diagrams for the different incoming states gg, qq, gq, qg are given in 



Fig. 4.6 d. The angular integrals have been calculated analytically, which leads to an analytic cancelation 
of the IR poles in e between the virtual correction, the real correction, and the mass factorization matrix 
elements squared. The latter one is described in section 2.2. % 

At lowest order, the cross sections for t\i\ and £2*2 production are given by the same analytical 
expression, since the mixing angle does not occur. At next-to-leading order the tig and four squark 
coupling introduce an explicit dependence on the mixing angle. The i^iz cross section can be obtained 



using the operation V\i described in eq.( |A.2| ). However, the dependence on the mixing angle turns out 
to be very weak. 

To perform a more detailed analysis the partonic cross section is expressed in form of scaling func- 
tions 



a 



'J 
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m~ 



/^(r ? )+47ra s (^ 2 



f7j +S {^ mj,e) + f$ (r,) + hM log I 



Ttl~ 



(4.5) 



where ij are the incoming partons, (77 = s/(4m~) — 1) with the partonic cm energy s, rrij genetically 
denote the set of masses entering the virtual corrections, and 9 is the stop mixing angle. For the sake of 
simplicity we have identified the renormalization and the factorization scale \it = fiR = \i. The scaling 
function f B contains the Born term, f v+s the virtual and soft-gluon contributions^ f H the hard-gluon 
contribution, and / the scale dependence. The function / combined with the running of the strong 
coupling constant and the scale dependence of the partem densities should yield a decreased dependence 
on /i. 

'Dividing the real gluon contribution into soft and hard gluons leads of course to some ambiguity in the definitions of f v+s 
and f H . This is in detail described in appendix |b| 
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The scaling functions are shown in Fig. \.\. The scaling functions f g q are identical to f qg . Only the 
scaling function f v+s depends on the mixing angle 9 and on the additional squark and gluino masses. 
The contribution of / + compared to f H is small in general, and the dependence on the mixing angle 
is suppressed. In contrast to the light-flavor squark production only the gluonic stop production cross 
section is proportional to 0, whereas the qq collision leading to an s channel gluon exhibits an over-all 

A Coulomb singularity similar to the one for the stop decays in 

0. 



factor 3 , as can be seen in eq.(4.2) 
section 



3TT| appears: The scaling function fY Q +S approaches a non-zero limit near threshold rj 



The emission of soft gluons from the incoming partons leads to an energy dependence log* near 
threshold. The leading log 2 terms are universal and could be exponentiated. All scaling functions 
approach a simple form in the limit -C 1: 
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(4.6) 



In the high energy limit rj 3> 1 the LO cross section scales ocl/s, eq.(42)- The NLO cross sections 
involving at least one gluon in the initial state approach a finite value in this limit. This is caused by the 
exchange of soft gluons in the t or u channel. Exploiting kx factorization [ i5 ] the non-zero limits of the 
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scaling functions can be determined: 
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The ratio of f gg and f qg is given by 2N : Cf which are the color factors for the exchange of a gluon 
between an incoming quark or one of the two incoming gluons, and the Born diagram, which in both 
cases is the gg one. 



4.1.3. Results 

The hadronic cross section is obtained from the partonic by convolution with the parton densities, 



eq.(2.5). The phase space integration for the Born cross section as well as for the real gluon emission is 



given in appendix B.l. 



Scale Dependence 

The dependence on the scale \x = fip = has been analyzed for the production of i\ pairs both at the 
upgraded Tevatron and at the LHC. The hadronic cross sections include a s and the CTEQ4 [|l7|] parton 
densities consistently in LO or NLO, which also enters the definition of the K factor K = (Tnlo/o'lo- 



The scale dependence is shown in Fig. 12. The leading order dependence is monotonic and varies by 



about ~ 100% for scales between mj/2 and 2mj .. The increase for small scales results from the large 
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Figure 4.2: Left: The renormalization/factorization scale dependence of the total cross section for t\ 
pair production at the Tevatron and the LHC. the maximum for the NLO cross section at the LHC is 
reached only for very small scales; Right: Effect of the variation of the scale on the upgraded Tevatron 
production cross section, as a function of the stop mass. The LO and NLO bands show the improvement 
of the theoretical uncertainties in the derivation of mass bounds. The SUSY scenario determining the 



masses is given in eq.( 1.27 ) 
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Figure 4.3: K factors for diagonal stop-pair production at the upgraded Tevatron and the LHC for a 



sample of stop masses. Scale choice: \i - 
initial-state gg and qq fractions are given. 



m 



tv 



For a comparison of Tevatron and LHC also the LO 



value of the running QCD coupling in leading order. At next-to-leading order the variation with the scale 
is reduced to ~ 30%. The monotonic behavior of the leading order curve is corrected in next-to-leading 



order, yielding a maximum value at some small scale, Fig. |4.2 



Supersymmetric Parameter Dependence 



The total hadronic cross sections at the Tevatron and at the LHC are given in Fig. 4.5. The masses in 



volved are fixed by a SUGRA inspired scenario eq.( |1.27| ). All squarks except the top squark are assumed 
to be mass degenerate. The mass range of the outgoing stop is varied independently of the other mass 
parameters. The same is done for the mixing angle and for the gluino mass. The dependence of the cross 
section on these internal mass parameters is small, as can be seen from the finite width of the central lines 
in the figures, i.e. the cross sections depend essentially on the outgoing masses and do not distinguish 
between t\ and £2- The light-flavor squark and the gluino contributions appearing in the loops are de- 
coupled even for numerically small masses. The search for stops therefore yields limits on their masses 
independent of any other parameter, unlike the squark/gluino or the neutralino/chargino case. 

One exception of this behavior is the kink in the next-to-leading order cross sections for the heavy 
stops. Similarly to the decay width of a heavy stop to a gluino, Fig. [33| , threshold contributions occur. 
In the stop production case the heavy stop can decay into an on-shell gluino and a top quark. The kink 
will be regularized by introducing a finite width for the stop, and for an analysis it has to be removed by 
resummation. However, for the search for stops at the Tevatron this parameter region is not of interest. 

The strong mass dependence of the K factor is due to different K factors for the quark and the gluon 
channel, both of which are only weakly mass dependent. However the contribution of the two channels 



varies strongly with the mass of the external particle, as can be read off Fig. 4.3. Whereas the gluonic K 
factor is large (about 1.3), the quark K factor tends to be smaller than and close to one. Weighted with 



the fraction of the incoming state these combine to the K factor given in Tab. |4.3[ Since at the Tevatron 
the contribution of the incoming quarks decreases as the stops become light, the K factor grows from 
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1 .03 to 1 .43 far from threshold. At the LHC the gluons dominate over the whole considered mass region 
and the K factor varies between 1.27 and 1.50. 

With cross sections between 0.1 and 100 pb the integrated luminosity f C = 20fb~ 1 should be 
sufficient for collecting a sample of 10 3 and 10 6 stop events, provided the particle exists and with a mass 
less than 450 GeV. The LHC with an integrated luminosity of / C = 300fb _1 could collect 10 5 to 10 8 
stop events in the mass range of 200 to 500 GeV. 

The normalized differential cross sections with respect to the transverse momentum and the rapidity 
is shown in Fig. [k4j. The transverse momentum of the outgoing stops is shifted to a softer regime by the 
momentum carried by the additional jet in the final state. A naive description using the K factor would 
not take into account this shift and therefore lead to large errors in the size of the cross section for a 
certain value of pr- The rapidity distribution keeps almost the same shape in NLO as in LO. However, it 
is not symmetric in NLO anymore. 



4.2. Non-diagonal Stop Production 



In hadron collisions mixed pairs t\t2 and t^i cannot be produced in lowest order, unlike in e + e~ colli- 
sions, since the involved coupling conserve the chirality eigenstate, which does not hold for the coupling 
to a Z. The mixed production cross section is therefore 0(af ). For a general mass scenario it is small 
but difficult to calculate. In the diagonal production we observe that the limit of a decoupled gluino gives 
a good approximation for the size of the cross sections. Therefore we calculate the production cross 
section for i\t2 + £2*1 in this limit. Only two one-loop diagrams contribute to the amplitude in the limit, 



they are given in Fig. fk6| e. They involve the production of diagonal stop pairs in gg fusion, which are 
rescattered to mixed pairs by the four squark vertex. The incoming quarks are suppressed. 
The evaluation of the loops yields 



Ooo[gg -> *i*2 + hh] = sin 2 (40) 
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13824 2vrs 3 



m~ log 2 (xi) 



m~ log 2 (x 2 ) 



(4.8) 



where the subscript in the cross section indicates the limit m- g — > oo. The coefficient A 1 / 2 is the usual 



2 -particle phase-space factor, i.e. A = [s — (m^ +mj ) ][s — (mj 
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CT [fb] 


_ limit 
Uqq V q q 


limit 

"39 "gg 


Tevatron 




0.201 -10 3 0.202 -10 3 
0.333 0.337 



0.087 -10 3 0.087 -10 3 
0.016 0.016 

0.131 -10" 4 


LHC 


t\t2 + titl 


4.137 -10 3 4.150 -10 3 
0.169 10 3 0.172 -10 3 



70.13 -10 3 75.00 -10 3 
1.422 -10 3 1.458 -10 3 
0.149 



Table 4. 1 : Cross sections for diagonal and non-diagonal pair production at the Tevatron and the LHC, 
using the default SUGRA-inspired scenario eq.(1.21). The non-diagonal results are given without the 
mixing factor sin 2 (A9).The superscript 'limit' denotes the asymptotic value of the cross section for large 
gluino masses. 
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the logarithmic discontinuities are defined properly by the infinitesimal shift s — > s + ie in fy. The 
fraction 37/13824 originates from the color factor (N - 1)[5N 2 - 2(N - l) 2 ] / [2567V 3 (N + 1)]. 

The cross section depends strongly on the mixing angle 8 through the overall factor sin 2 (46*). Nu- 
merical values for the diagonal and non-diagonal pair cross sections are compared in Tab. [4.1[ Note that 
the mixed-pair cross section is given in this table without the mixing factor sin 2 (40). The values for the 
cross section for producing mixed stop pairs in the large m§ limit are very small at the Tevatron as well 
as at the LHC. 
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Figure 4.4: The normalized differential cross section for the production oft\ pairs at the Tevatron. The 
mass scenario using the central scale is defined in eq.( 1.2^ ). 
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Figure 4.5: The total cross section for t\ and t<i pair production at the Tevatron and at the LHC. The band 
for the NLO results indicates the uncertainty due to the scale dependence. The mass scenario is given in 



eq.< \1.2'fy . The line thickness of the NLO curves represents the variation of the gluino mass between 280 
and 900 GeV and of the mixing angle over its full range. The kink in both of the cross section results 
from the on-shell decay of the heavy stop to a very light gluino and a top, and can be regularized by 
introducing a finite stop width. 



44 



CHAPTER 4. PRODUCTION OF SCALAR TOP QUARKS 




(c) 



etc. 




Figure 4.6: Generic Feynman diagrams for the diagonal stop pair production: (a) Born diagrams for 
quark and gluon incoming state; (b) vertex corrections; (c) box contributions; (d) real gluon/quark 
emission for different incoming states; (e) mixed stop pair production in the limit of decoupled gluinos. 
The self energy contributions are not shown. 



5. R Parity Violating S quarks 



5.1. Production in ep Collisions 



Limits on the mass and the coupling A' as defined in eq.(1.6) can be derived from the direct search at 
different colliders. HERA ep scattering could produce squarks via the R parity violating coupling A' to 
quarks and electrons, where the flavor of the squark has to be chosen consistent with the whole set of 



current bounds [ 53 , 24 , 25 ] : 



eq 



(5.1) 



This resonant s channel production process can be described in terms of general scalar leptoquarks. The 
Yukawa matrix A' does not have to be diagonal in flavor or generations. 

For the leading order hadronic cross section the convolution with the parton density, as defined in 
eq.(^5) becomes trivial, since the energy-momentum conservation yields a factor 5(1 — m 2 /(xS)) 



i2 



0"LO 



4m 2 

2 



m 



J a 



m 
~S 



M C"LO 



(5.2) 



The parton densities are taken at the factorization scale \i, and possible flavors of incoming quark are 
fixed by charge conservation and the charge of the outgoing squark. 

The NLO contributions consist of virtual gluons, real gluon emission, and the crossed eg incoming 



state. All other supersymmetric particles in this scenario are assumed to be decoupled, see appendix A. 3. 



Some generic Feynman diagrams for the matrix element can be derived from Fig. [^J by replacing the 
external gluino by a positron and removing all internal gluino contributions^. After renormalization and 
mass factorization the dimensionally regularized NLO cross section is finite and can be written as 



m ePf \ * 
f NLO = ~g Jq (X) OLO 

A 



Cfoi s 
-i 42 



7T 



+ A + A c 



^ a LO dy /f (y) j - Z -P qq {z) log z + C F (1 + z) 



+ C F 



' log(l - z) 
1 - z 



l-z 



(2 + z + z 2 )log(l-z) 



A g = ^a LO fdyf q p (y)d-P q9 (z) 



log- 



l-z) 



+ 2 



+ z(l — z) log z + 1 



(5.3) 



'The next-to-leading order calculation was performed in parallel to [p^l, and the results are in agreement 
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where x = m 2 /S and z = x/y, and the renormalization and factorization scales are set /xp = fiR 



m. 



The splitting functions Pij are defined in eq.( B.23 ), and the + distributions are given in eq.( B.24 ). The 
NLO cross section is proportional to <3"lo> i- e - tne anomalous coupling drops out of the K factor. In the 
formulae the factorization and renormalization scales have been identified with the squark mass. 

The 



The K factors for the production of an up-type and a down-type squark are given in Fig. |5T 



next-to-leading order correction for the HERA process is dominated by the gluon emission from an 
incoming quark. They are positive in the whole mass region considered. The virtual correction as well 
as the crossed channel are suppressed. Since the parton distributions enter the K factor, the different up 
and down type valence and sea quarks receive different corrections. 



The dependence on the factorization and renormalization scale can be made explicit in eq.(5.3) by 
adding log n 2 FR /m 2 to the Q2 term in the virtual correction and to the logarithms multiplied by P^. 
Additionally the mass m in the argument of the parton densities and the running coupling A' has to be 



replaced by fj,. The decrease of the scale dependence of the cross section is shown in Fig. 5.1 



The partial width for the decay of a R parity violating squark can be computed using the results for 
the scalar top squark. The width for the decay channel 



eq 



(5.4) 



in NLO reads 



A m 
16vr 



1 + 



27 



2C2 



CFOi s 



7T 



(5.5) 



This correction is small (~ 10%). If the R parity breaking coupling A' is small, then the squark will be 
long-lived, and the NLO correction will not change this effect. A typical scenario could be a squark with 
a mass of 200 GeV and a Yukawa coupling A' ~ e/10, which leads to a decay width of V ~ 3 GeV. 
Bound states however will only occur, if R parity conserving decay modes are kinematically forbidden. 



47 



5.2. Production in Hadron Collisions 

In contrast to the HERA production process the production of squarks at hadron colliders, R parity 
conserving and R parity violating, is fixed by the QCD coupling, as long as the production process does 
not involve any weakly interacting particles^ In the scenario under consideration all other non-Standard 
Model particles are decoupled. The Feynman diagrams for the Born cross section [|l4|] are the same as 



for the stop pair production, Fig. 4.6a 



ml 99 



(5.6) 



The Born cross section is given by eq.(|4.2|), the t channel diagram for quark-antiquark collisions is 
in this case not suppressed by the incoming state but absent, due to the decoupling of the gluino. The 
strong coupling is independent of the flavor of the light squark, i.e. in case of hadroproduction all flavors 



look identical, unlike Fig. 5.1. The calculation for R parity breaking squark production is the same as 
for the stop pairs, the Feynman diagrams are given in Fig. ^ taking the limit of decoupling gluino 
and removing the strong four-squark coupling. This follows from appendix [A|. The partonic NLO cross 
section only depends on the partonic cm energy s and on the mass of the squark, leading to the [in this 
case literal] scaling functions for the Born, virtual and soft, hard, and scale dependent contributions 



or 



^(r/)+47ra s (// 2 



v+s 

ij 



(5.7) 



where i, j denote the initial state partons. The factorization and the renormalization scale have been 
identified, and rj = s/(4m~) — 1. The numerical form and the structure of these scaling functions are 
very similar to the stop case, Fig. |4.l| , since the SUSY parameters except for the outgoing mass influence 
the stop scaling functions only marginally, and the four-squark vertex only contributes to the numerically 
suppressed virtual corrections. This result for the R parity violating squarks can also be obtained from 
the even more general light-flavor squark production M\ in the limit of a large gluino mass. However, 



2 Other production channels [ p7| involve the weak coupling constant and more than two final state particles; their cross 
sections are significantly smaller than the squark-antisquark production via an s channel gluon. 
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Figure 5.2: Renormalization/factorization scale dependence of the total cross section pp — > qq at the 
non-upgraded Tevatron. The arrow indicated the average invariant energy < s y 1 / 2 in the hard subpro- 
cess, which was used in the original analysis. 
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Figure 5.3: Total hadronic cross section for the production of R parity violating squarks in hadron 
collisions pp — > qq. The leading order result is given for the renormalization factorization scale p = m. 
In the first case also the values of a s and the parton densities have been taken in leading order. The K 
factor for the smaller central scale m is comparably small. 



the occurrence of the gluino mass in the Born term requires a numerically large gluino mass for the 
analysis. And the four-squark coupling has to be removed as for the derivation from the diagonal stop 
pair production. 

Exactly as for the stops, the scale dependence Fig. 5^ in NLO leads to a maximum and an increasing 
accuracy for the derivation of limits on the mass of the particles from non-observation. Since in leading 
order any scale of the process can be considered, choosing the invariant energy of the final state p = yfs 
is possible. Especially if NLO calculations are not available the choice of the scale of the process leads 
to considerable uncertainties, which is illustrated in this example: In NLO the choice of the factorization 
scale is no longer free, since it must be a parameter defined in terms of external variables, which does 
not allow for p = yfs. The choice of p = yfs in next-to-leading order leads to an inconsistence of the 
order a s , independent of the order of perturbation theory, in which the process has been calculated [48]. 
Another scale for the production of massive particles would be the final state mass m, which turns out 
to be much smaller than y/s averaged with the weight of the cross section. The different K factors for 
both choices of the scale are given in Fig. 5_2. For the typical behavior of the leading and next-to-leading 
order cross section the large scale, y/s, is far from the point of maximal convergence. This reflects the 
appearance of logarithms log m 2 /p 2 , which render the corrections unnaturally large. 

The NLO cross section and the LO cross section for the central mass scale for the upgraded Tevatron 
are given in Fig. 5.3. The variation of the cross section with the scale decreases from 100% to 30% if 
the NLO result is used. Compared to the p = m leading order result this leads to a more accurate and 
always higher mass limit, the improvement of the mass bounds however stays below 10 GeV for the 
central scale. 



Experimental Analyses 

One of the design features of the ep collider HERA is the search for leptoquark-like particles, i.e. particles 
which carry electron and quark quantum numbers, and can in case of scalars be identified with R parity 
violating squarks. They can occur either in the s or in the t channel and lead to an excess in the ep cross 
section [49]. The decay channel of this squark strongly depends on the masses of the supersymmetric 
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scenario considered, but may well include a high-p^ electron. This signal is essentially background 



free [|23|]. The interpretation of the combined data of ZEUS and HI as R parity violating squarks leads to 
combined limits on the Yukawa coupling A' and the branching ratio BR to the observed high-p^ electron 
final state [The numbers are based on the data analyzed by fall 1997.]. 
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(5.8) 



Couplings to d and u would lead to a large enhancement in the e~p run and are forbidden because of 
their non-observation. Moreover, couplings involving a positron and an up-type quarks would lead to 
electric charges, which do not occur for MSSM-type squarks. The coupling A' has to be interpreted 
as entries into the non-diagonal Yukawa coupling matrix, connecting down-type quarks to squarks and 
electrons. The diagonal matrix element A' m , which would lead to the production e + d — > u, is excluded 
by neutrinoless double beta decay. Possible candidates for a resonance production are 



~d -> c L ,t L 



e s 



t L 



(5.9) 



As depicted in section |1.2.1| , atomic parity violation yields strong limits on X' 2 /m 2 for any lep- 
toquark interacting with valence quarks. They can for a mass of ~ 200 GeV be translated into 



bounds on the Yukawa coupling matrix A' < 0.055 [e + d] [ )24| , P5Q . Combined with the measured val- 
ues at HERA, eq.(|j|), this yields lower limits on the branching ratio to the observed eq final state 
BR > 0.2 • • • 0.4 [e + d]. The limits obtained from atomic parity violation are derived for the presence of 
only one particle being responsible for the possible deviation from the Standard Model. More than one 
R parity violating squark influences this analysis, the result depending on the sign of interference terms 
and thereby on the quantum numbers. The assumed left handed stop quark is in general a superposition 
of two states with different mass and equal electroweak properties. This strengthens the bound on the 
branching ratio: 

/ 2 \ 

/ m ^ \ 



BR — > BR ^ 1 + tan 2 9 — Jr ) (5.10) 



At LEP, the obtained limits on A' are relevant for sea quarks only [24, E5h. By the same token as for 



atomic parity violation they start from A' < 0.6 [e + s] and give BR > 0.05 • • • 0.2 [e + s]. 

The data from the search at the Tevatron can be written — given the mass of ~ 200 GeV from the 
HERA analysis — as an upper bound in the branching ratio BR < 0.5 • • • 0.7. Theoretically the compe- 
tition between supersymmetric R parity conserving, and R parity violating decays makes it possible to 
vary the branching ratio into the eq mode with the mass of the particles forming the decay chains. 

In the case of ed — > cl the most important MSSM-like decay modes are cx° and s\j ■ Assuming the 

gauginos being heavy [m ~+ > 200 GeV] insures that the branching ratio lies in the region of BR ~ 1/2. 

1 

The higgsino decays suffer from the small strange quark mass. For ed — > II the strong decay channels 
are forbidden, as depicted in section ||. The decays into tx® and bxlj can be suppressed by large gaugino 



and higgsino mass parameters. However, suppressing the whole set of possible stop decays, eq.( |3.1[ ), 
yields some fine tuning of masses and mixing also in the sbottom sector. The coupling A' for the es 
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tL production channel is comparably large, which renders the different Yukawa coupling and gauge 
coupling mediated decay channels of similar size, and thereby prevents from any fine tuning. 

The interpretation of the HERA excess as R parity violating squarks is therefore not ruled out by the 
bounds set by other experimental and theoretical analyses, but give an impression how different collider 
experiments and non-collider experiments like atomic parity violation and search for neutrinoless beta 
decay can altogether constrain the parameters in the same model. Incorporating all available data, HERA 
itself is the only experiment able to remove this interpretation by non-confirming the excess. 



6. Conclusions 



In this work supersymmetric QCD corrections to decays involving scalar top quarks and to the hadropro- 
duction of neutralinos/charginos and stops are presented. The decay widths as well as the production 
cross sections calculated in perturbation theory exhibit an unphysical dependence on the renormalization 
and/or factorization scales. In leading order this dependence is in general strong. Compared to a central 
scale, which could be the mass m for the decay width and the hadro-production cross section of massive 
particles, this leads to variations up to a factor of two for scales between m/2 and 2m. In next-to-leading 
order this dependence is considerably weaker, i.e. in addition to the K factor, the next-to-leading order 
results always improve the precision of the theoretical prediction used for the experimental analysis. 

The basic properties of the scalar top sector are investigated by calculating the supersymmetric QCD 
corrections to the strongly and weakly coupling decays, including a stop either in the initial or in the final 
state. An elegant definition of a running mixing angle in next-to-leading order is given, in order to restore 
the Born type symmetries between the stops in next-to-leading order observables. The mixing angle 
counter term is compared to other renormalization schemes. Although the phenomenological motivations 
for the various schemes are different, the numerical differences are shown to be small. 

The different stop decay widths obey a strong hierarchy, starting with rare decay channels, and then 
proceeding towards weak and strong two-body decays for an increasing mass of the decaying stop. The 
strong decay will be dominant for a heavy stop state. The next-to-leading order corrections to the heavy 
stop state decaying into a gluino are large ~ 30% and always positive, while for the gluino decay into 
a top squark they turn out to be small and negative ~ —5%. This feature also arises for the light- 
flavor squarks and is due to interference between different color structures and the different analytical 
continuation of logarithms. The dependence of all decay widths on the mixing angle can be described by 
a K factor, which stays constant for varying angle. The scale dependence of the decay widths is reduced 
from a factor of two to about 50% in next-to-leading order. 

The weak decays of a light stop into a neutralino and a chargino are analyzed, to illustrate the running 
mixing angle. The next-to-leading order supersymmetric QCD corrections are small compared with 
typical strong decays. Apart from special mass scenarios and threshold effects they are < 10%. In 
contrast to the strong decays the sign of the correction is not fixed, it is strongly dependent on the mass 
scenario considered. The same holds for decays of heavy neutralinos. They can be produced at e + e~ 
linear colliders and will in supergravity inspired scenarios be higgsino-like. Therefore the decay induced 
by the top-stop Yukawa coupling can give large contributions, whereas the light-flavor final states are 
strongly suppressed. The next-to-leading order corrections to these widths are moderate: ~ 10%. 

The next-to-leading order production cross section for neutralinos and charginos can be used to derive 
mass limits at the upgraded Tevatron and at the LHC. This yields an improvement of the mass bounds for 
these particles obtained at LEP2. Although mass and mixing parameters could be derived from cascade 
decays of strongly interacting particles, the only way to keep maximal independence of the choice of the 
model is the direct search. Similar to the case of gluino production, on-shell and off-shell intermediate 
particle contributions have to be distinguished. This is done in a manner, which naturally coincides with 



52 



CHAPTER 6. CONCLUSIONS 



the experimental analyses. The different higgsino-like and gaugino-like contributions to the production 
cross section can be analyzed and give a smooth picture of the next-to-leading order corrections. Al- 
though the scale dependence in next-to-leading order is not as much improved as for strongly interacting 
particles, where the dependence on the running QCD coupling arises in leading order, it stays below few 
percent in next-to-leading order. The K factor for all possible final state neutralinos/charginos is nearly 
constant for varying masses and corrects the leading order result by +20% to +50%. However, strong 
cancelations between different diagrams may lead to large K factors in the mixed production channel, 
strongly dependent on the mass and mixing parameters chosen. 

The search for scalar top quarks is naturally the next step after the search for light-flavor squarks 
and gluinos at the upgraded Tevatron and at the LHC. Since the QCD type couplings are invariant under 
chiral transformations, and cannot distinguish between the right and the left stop, the production cross 
section in leading order depends only on the mass of the produced particles. The mixing angle as well as 
the mass of the light-flavor squarks and the gluino only enters through the virtual corrections. This de- 
pendence is found to be much smaller than the scale dependence and thereby negligible. The corrections 
to the diagonal stop pair production are different for quark and gluon incoming states. Whereas they are 
small and negative ~ —5% for incoming quarks, they are large and positive ~ 50% for gluons. Since 
even at the Tevatron only very heavy stops are produced mainly in quark-antiquark collisions, this mass 
dependent K factor leads to an increase of the mass bounds derived in leading order. In addition to the 
total cross section the differential cross sections are analyzed: The rapidity distribution at the upgraded 
Tevatron can be described using a K factor, the transverse momentum distribution however is shifted to 
a softer regime in next-to-leading order. 

The search for mixed stop final states could be a means to measure the stop mixing angle directly 
at hadron colliders. However, the cross section is based on a one-loop amplitude, and the calculation in 
the limit of decoupled gluinos shows, that it is much smaller than the cross section for the production of 
two heavy stop particles. The direct measurement of the stop mixing angle will be possible only in the 
process e + e~ — * tit>z. 

The calculation of the stop production cross section can be adapted to the search for supersymmetric 
R parity violating squarks at the Tevatron. Similar to the HERA production process, the next-to-leading 
order calculation for the hadroproduction gives rise to improved bounds on the mass and on the branching 
ratio into the observed eq channel. However, this is not sufficient to close the window in the branching 
ratio between the atomic parity violation and LEP on one side and the Tevatron on the other side, before 
the excess completely vanishes by accumulating more HERA data. 

All calculations have been performed in a supergravity inspired GUT scenario. This connects not only 
the different gaugino masses, but also the masses of the squarks and the stops. All these particles can be 
searched for at the upgraded Tevatron and at the LHC. In Fig. ^[l| the relevant cross sections are given as 
a function of the mass of the produced particles. The weakly interacting particles are strongly suppressed 
at the LHC; however, the search for leptonic events at hadron colliders is completely different from the 
hadronic final states. Assuming a similar efficiency for the search for stops and light-flavor squarks the 
search for both of them in parallel seems to be promising, in particular, since the search for squarks 
is a two dimensional, the search for stops a one dimensional problem. The smaller stop cross section, 
due to the missing t channel gluino contribution and the fixed non-degenerate flavor, is compensated 
by the small mass of the light stop state in typical scenarios. For gluinos, like-sign leptons in the final 
state should improve the efficiency considerably, yielding all production processes given in Fig. very 
promising at the upgraded Tevatron and the LHC. 
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Figure 6. 1 : The total cross section for pairs of squarks, gluinos, stops, and neutralinos/charginos 
as a function of the mass of the produced particles. The usual SUGRA mass spectrum mo = 
100 GeV, m-i/2 = 150 GeV, Aq = 300 GeV, tan (3 = 4, fi > is denoted by the arrows. The masses of 
X2 and X\ differ only by a few GeV in the small regime. For the neutralino/chargino production 
cross section the masses are consistently varied with m 1( / 2 . The cross sections are given for the upgraded 
Tevatron and for the LHC. 
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A. SUSY Lagrangean 



A.l. Feynman Rules for Supersymmetric QCD 



In this appendix we give a complete set of Feynman rules, as used in the calculation of the various 
processes. Since fermion number violating processes have to be considered, the rules make use of a 



continuous fermion flow [50], which has to be fixed once for any process. The Dirac trace has to be 
evaluated in the opposite direction of this fermion flow. If not stated otherwise, we assume the fermion 
flow being identical to the Dirac fermion flow e.g. for the quarks. 

Using these Feynman rules makes it needless to introduce charge conjugation matrices. Moreover, the 
relative signs of different diagrams contributing to the same process can be fixed: Any permutation of 
the fermion flow of two external fermion lines gives rise to a factor (—1) for the matrix element, due to 
Fermi statistics and Wick ordering. To match the spinors for processes including two external Majorana 
and two external Dirac fermions, all diagrams have to be evaluated with two different directions of the 
fermion flow. We have checked explicitly that the two possible orientations of the fermion flow of two 
combined diagrams lead to the same result. 

By changing the orientation of the fermion flow, the signs of vertices with different types of couplings 
[S, P, V, A] change: 



coupling 




sign 


scalar (5) 


cic- 1 = 


1 T 
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pseudo-scalar (P) 


c^c- 1 = 


T 

7 5 
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vector (V) 


c^c- 1 = 


T 




axial-vector (A) 


C (757m) C- 1 = 


(757/J T 


+ 



Table A.l: Transformation of couplings with the orientation of the fermion flow. C is the charge conju- 
gation matrix. 



Standard Model Feynman Rules 

All momenta in the Feynman rules are defined incoming. The Standard Model couplings of quarks, 
gluons, Fadeev-Popov ghosts and weak gauge bosons are given in Fig. Al in the Feynman gauge. The 
generators of SU(3)c obey the relations 



Tr T a T 



-5 ab 
2 



rpa rj-it 



if ahc T c 



(A.l) 



The fermion propagators are defined as i/(jf> — m + ie), where p is the momentum in the direction of the 
fermion flow. The fermion number flow does not occur. The gluon propagator is —ig^j (p 2 + ie). 



56 



APPENDIX A. SUSY LAGRANGEAN 




Figure A. 1 : Feynman rules for Standard Model quarks, gluons, ghosts and weak gauge bosons, the 



dotted lines for scalars refer to Fadeev-Popov ghosts. The generic couplings £, r are defined in Tab. 4.2 
Plr a re the chiral projectors (1 =p 75) /2 



Supersymmetric QCD Feynman Rules 



The Feynman rules for supersymmetric QCD include, besides the Standard Model particles, the gluino 
(g), the light flavor squarks (q L , q R ), and the mixing top squarks (ti, £2)- Like in Fig. |A.lj all momenta 
are defined incoming. The coupling qqg preserves the helicity of the quark and its scalar partner as 
well as the flavor. In higher orders it has to be modified to restore the supersymmetric Ward identity, as 



described in chapter |L5J. The same holds for the qqg coupling. In leading order we use g s for all strong 
Standard Model and their supersymmetry transformed couplings. 

The qqg(g) vertices preserve the flavor and the 'helicity' of the squark. Since only two squarks are 
present, this coupling cannot mix the mass eigenstates. 

The qqg vertex Feynman rule is given for a light-flavor q LR with the corresponding sign and projector. 
To obtain the rules for i\ one has to add the contributions of the helicity eigenstates and multiply the L 
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Table A. 2: Couplings of quarks to weak gauge bosons as used in Fig. |4.7| . 
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Figure A. 2: Feynman rules for supersymmetric QCD. The dotted lines denote the squark q L ,q~R- The 



tensors are defined in Tab. \.3; Plr are the chiral projectors (1 =p 7s)/2 



term with sin 8 and the R term with cos 6. The mixing of the scalar top quark is described in detail in 



chapter 1.2.4 The relevant terms for stop mixing in the Lagrangean include the coupling to the gluino 
and the four squark coupling which arises from the D term in the scalar potential, described in eq.(L4)]. 
They can be expressed using the permutation operator 



Vn ■ [h *-> h; Cq -> -s§, Sg -> c- B 



(A.2) 



which links the vertices including t\ and t% 



V2g s T% (1 + Vi 2 )h 



cz Pl ~ sz Pr 



tj h* + h.c. 



1 + Vvi) ti* hj L 2 2§ Sf l hi in + 2 [s 2 2§ Sf l - Sf l 



+ 4 c 2§ S\ lkl {qlI qa - q R * k q R i) 



(A3) 



The sine/cosine of an angle (3 is as usually denoted by sp, cp. However, these vertices describe only 
processes which are essentially diagonal either in i\ or in t^. Four squark vertices, mixing i\ and I2, can 



'This coupling is not fixed by the requirement that the scalar top quark should carry a fundamental representation SU(3) 
charge. 
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be derived from the Lagrangean 

,2 



r' 



9, 



4 S 29 hj + *2* hj) |c 2 e ^ *li _ h*k hi) 

+ 2 s{ jkl (&* fe; - fefc g«z) | 



(A.4) 



The structure of the four squark coupling is given in terms of the flavor fj and the helicity of the Standard 
Model partner hj. The two tensors used in Tab. A3 



are 



S 



S 



(i) 

ijkl 

(2) 
ijkl 



5il$jk — JjSijfikl 



N - 1 
N 



(Sirfjk + Sij5 k i) 



(A.5) 



mjMi 


Sijkl 


qqqq [hi = 


hj hj = hi fi = fj fj = fi] 


o(2) 
D ijkl 


qqqq [hi = 


hj hj = h t fi + fj fj / fi] 


<j(l) 
°ijkl 


qqqq [hi / hj hj ^ hi] 


^°ijkl 


hhQLRQLR 




~r c 29 D ijkl 


hh^LR^LR 




±c - S W 
^29 °ijkl 


hhQLRQLR 


hhQLRQLR 


±s - S (1) 

3:6 20 °ijkl 


hhhh 


hhhh 


c 2 9 (2) 


hhhh 


hhhh 


+s -c - S {2) 

^ b 29 c 29 °ijkl 


hhhh 


hhhh 


-s -c - S {2) 

b 29 c 29 °ijkl 


hhhh 




s 2 c-(2) 


hhhh 


hhhh 


,2 q(2) 



Table A. 3: Tensors arising in the generic four squark coupling in Fig. A.2. The tensors and are 
defined in eq. f |4.5| ). 



A.2. Neutralinos and Charginos 
Diagonalization for Neutralinos 



The diagonalization procedure for Neutralinos is described e.g. in [plQ. The four Majorana neutralinos 
(Xj 1 ) are defined as the mass basis of all neutral higgsinos, the photino, and the zino interaction eigen- 
states. They are four component Majorana spinors, therefore the mass matrix is symmetric. It is possible 
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to start from a non-diagonal matrix in the (BW3) basis or in the (jZ) basis. We denote the two possible 
mass matrices as M and M.' for the two component states: 



M 





-m z s w C/3 
\ m z s w Sj3 





m z c w C/3 
-m z c w sp 



-m z s w cp 
m z c w C/3 




m z s w Sf3 \ 
-m z c w Sj3 

/ 



(A.6) 



The corresponding unitary mixing matrices N, N' diagonalize these mass matrices, N in the (BW3) and 
N' in the (•yZ) basisQ. In case the mass matrix is real, the mixing matrices N, N' are also chosen to 
be real, to keep the couplings from becoming complex i.e. in this special case the in general complex 
unitary transformation becomes real and orthogonal. However, for CP invariant observables the typical 
coupling factors must be purely real or imaginary for any complex unitary mixing matrix. 



AT* M N- 1 = M diag 
N'* M' N'- 1 = M 



diag 



N'ji = Njic w + N j2 s w 
N'32 = -N jlSw + N j2 c v 
N' 



N 1 



J'3 
J'4 



^'3 

N j4 



(A.7) 



2 Any complex symmetric matrix can be diagonalized to a real diagonal matrix by a unitary transformation U T AU where 
= I/ -1 . The diagonal matrix can be chosen real since phase factors can be absorbed into the unitary mixing matrix. 




-t 7 " (LP L + RP R ) 



(LR) 




iV2(±Al R P LR + B* LR P RL ) 



(LR) 




iV2(±A LR P RL + B LR P LR ) 



Figure A. 3: Feynman rules for neutralinos and charginos. The dotted lines denote a left or right scalar 
quarks q L , q R . L, R are defined in Tab. AA A, B in Tab. |4.5| Note that the qqx vertices for incoming 
and outgoing squarks are linked. 
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Xi Xj 1 



xtxj z 



x°ixlw; 



xU~w+ 



L 



of i3 of * 3 - or^of * 



e5ij 



S7/.C7, 



V iX V 3 \ + -V l2 V* 2 - 5 l3 sl 
V2s w V 

(of * 4 ^- 2 - ( av> rt* a + Cm ^; 2 ) ^1 



R 

L* [Of fc4 <-> -Of fc3 ] 

L* [V -> E7] 
L* [F 17, Of i4 -» -Of 
L* [y -> 17, Of i4 -> -Of 



Table A. 4: Couplings of neutralinos and charginos to weak gauge bosons as used in Fig. \.3 The mixing 
matrix Of is defined in the photino-zino basis. For charginos, the fermion flow is assumed to follow the 
xf fermion number flow. 



The mixing matrix is defined in terms of the arbitrary-sign eigenvalues of the mass matrix, m, : 



Cm Tflf 



mi 



Nj2^ 

Nn s w rriw - mi 

N i3 _ (j,(m B - mi)(m^ - m») - m\spcp\{m B - m^)c 2 w + - rm] 
Ni ~ 



N4 

Ni 



m z s w {mw - mi)(p,cp + rmsp) 
mi{m B - mi){mw - m) - m 2 z c 2 p[(m B - m^)<? w + - mi] 



m z s w (mw ~ mi){ncp + rrusp) 



N, 



il 



1 + 



N2 
Ni 



+ 



N3 

Ni 



+ 



A^4 

Nu 



-1/2 



(A.8) 



The entries of .Mdiag are not necessarily positive, if the mixing matrix is kept real, i.e. the eigenvalues 
are only equal to the physical masses up to a sign. It is possible to work with a Lagrangean including 
negative mass eigenvalues rrij . In the final expression these eigenvalues have to be substituted by their 
absolute values rrij — > ±m^ , in order to express the analytical result in terms of physical masses. An 
equivalent way of introducing these phases is to define a complex mixing matrix [01, Of], in which a row 
is multiplied by i, if the corresponding eigenvalue is negative: 



N 



ki 



OTfc/ [I = 1, 4] if eigenvalue m& positive 
rJXki [I = 1, 4] if eigenvalue m& negative 



(A.9) 



The re-definition of N' — > Of is defined in analogy; in typical scenarios one of the higgsino eigenvalues 
rrik [k = 3, 4] turns out to be negative, whereas the re-rotation 0T — ► Of only affects the gaugino part of 
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the mass matrix. Using this matrix 91, 91' one can always stick to the positive mass values. 

m* M 9T 1 = Diag [m %j ] 
9f* M' = Diag (m^) j = 1, 2, 3, 4 (A.10) 

The masses of the four neutralinos are re-ordered by their size after diagonalization, where \\ is defined 
being the lightest of the four. Combining the complex couplings including the matrix 91, 91' leads to 
exactly the same analytical results for CP invariant observables as using the matrix N, N'; the phase 
factors from the negative masses now enter by collecting factors of i 2 in the typical combinations of the 
couplings and by anti-commuting the Dirac matrices. One advantage of the latter strategy is, that the 
neutralino mass matrix is not fixed to real values by first principles [|52|]. 



Diagonalization for Charginos 

Charginos (xJiXj) are the mass eigenstates of charged winos and higgsinos. The positive and negative 
charge particles mix independently. Since the charginos are no Majorana particles the mass matrix is not 
symmetric. Nevertheless the Dirac-chargino vertices can be fermion number violating. 

% V2m w S{3 \ (All) 

The unitary diagonalization matrices for the positive and negative winos and higgsinos are V and U, and 
the eigenvalues of the diagonalized mass matrix can in general assumed to be real[| The mixing matrices 
themselves are only real, if \i is chosen to be real. 

U* M V' 1 = Diag (mfc) 3 = 1,2 (A.12) 




Neutralino/Chargino Feynman Rules 



The Feynman rules for the neutralinos and charginos are given in Fig. A3. The qqx vertex is given for 
left and right squarks, the coupling to the mixing scalar top quark is a superposition of both couplings, 
as it is for the gluino case. 



A.3. R Parity breaking Squarks 



The breaking of R parity only adds new interaction terms to the Lagrangean, eq.(1.6), which are not 
related to any of the Standard Model gauge symmetries. For any of these scenarios including non-MSSM 
squark couplings, the QCD Feynman rules are still completely fixed by the requirement, that any squark 



should be part of the fundamental representation of SU(3)c, i.e. carries quark-type color charge Q53Q. 
Under the simplifying assumption of one light squark flavor, one can integrate out all heavy strongly 
interacting supersymmetric particles and regard the R parity violating scenario as the extension of the 
Standard Model by one leptoquark-like squark. Another difference between a scalar leptoquark model 
and the MSSM squark sector occurs: The MSSM four-squark coupling originates from the D terms in 



the scalar potential eq.(1.4). In the most general effective model this term need not be present. The 



3 Any complex matrix can be diagonalized to a real and positive diagonal matrix using two unitary matrices U T AV. If the 
matrix A is real the matrices U and V can be chosen real. 
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Ar 




Br 


uux°j 

ddx° 

du (xj) 
\ J / out 

ud (xt) 

V J / out 


gs w Q<rtji + —m' j2 (T 3 - q4) 

9 TJ 


Al [T 3 = 0] 
A* L [T 3 = 0] 






2m w S/3 j4 

9 m d m ,* 

2m w Cj3 j 

2m w Sj3 32 

gm d LT * 

2m w Cf3 j2 


Bl 
Bl 

2m w cp 

gm u 

2m w S[3 



Table A. 5: Couplings of neutralinos and charginos to squarks and quarks as used in Fig. [4.j[ The mixing 
matrix 9t' is defined in the photino-zino basis.. 



four-squark coupling will then be proportional to the weak coupling constant, and quadratic divergences 
of scalar masses occur. 

The coupling to any pair of Dirac and Majorana fermions can be identified with the most general param- 
eterization of the t\tx°j coupling. However, the qqe coupling constant is not fixed by any gauge coupling, 



but results from the superpotential term eq.(1.6) 



C int = X'eqq + h.c. (A. 13) 

Since its Dirac structure is fixed by the helicity of the quark and the electron, the coupling can be param- 
eterized as — i{AlPi + ArPr). For unpolarized particles this yields {A 2 L + A 2 R ) in the HERA production 
cross section and the decay width, likewise for leading and next-to-leading order, this means the coupling 
completely drops out of the next-to-leading order QCD K factors. 

Because the coupling A' connects two strongly interacting particles, it has to be renormalized, e.g. in 
MS. Furthermore, it runs in complete analogy to the running QCD coupling constant 

fM-^n (A , 4) 



B. Radiative Corrections 



B.l. Phase Space and Partonic Cross Sections 

(2 — >• 1) Production Cross Section 

Assume the production of one final state particle in NLO 

q(h) + q(k 2 ) -» Xtip^l+gih)] (B.l) 

with massless partons kj [j = 1,2,3] and a massive particle in the final state p\ = m 2 . The invariants are 
the usual Mandelstam variables 

s = 2{hk 2 ) h = 2(hk 3 ) u x = 2(k 2 k 3 ) (B.2) 

These Mandelstam variables can be expressed in terms of the rescaled gluon emission angle 8 

y = \(l + c e )e[0,l] 
t= - s(l - r)y 

u= -a(l-r)(l-y) (B.3) 



where r = m 2 / s and cq = cos 9. The formula for the cross section becomes []54[]: 



Kij are the spin-color averaging factors, for the HERA process K eq = l/(4iV). For the general squark 
production at HERA the integration over the gluon emission angle can be performed analytically. Poles 
in e appearing in the soft/collinear phase space region cancel with the divergences arising from the 
virtual gluon contributions and with the mass factorization, i.e. the renormalization of the parton densities 



described in section 2.2.2 



(1^2) Decay Widths 

The calculation of the different decays has been performed in four dimensions. The general decay mode 



is 



*i(pi) - X 2 (p 2 ) + X 3 ( P3 )[+g(k)] (B.5) 
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where all particles are assumed massive p 2 = m 2 [j = 1, 2, 3] and k 2 = A 2 . The Born and real gluon 
emission decay width is [] 



B 



t b = K 1 A^K ml^I) y 
16nm\ m\ 

rR ~ *>*k S § § S><- +- +- + *>E l""f 

iv"g = 1/iVis the spin-color averaging factor for the decay of the squark. The integrals 

= ~ 2 J 6 iPl+P2 +P ' + k) 2(k Pn )...2(kp Jm ) (B - 7) 

can be found in the literature [|55|]p|. 

(2^2) Production Cross Section 

The partonic production process of two massive particles can be written as 

q(ki) + q(k 2 ) - Xi(pi) + X 2 (p 2 )[+g(k 3 )} (B.8) 

The fcj [j = 1, 2, 3] are assumed massless, the pj [j = 1, 2] will in general have different masses m\,m 2 . 
Replacing quarks by gluons and vice versa does not have any effect on the kinematics. We introduce the 
invariants: 



s 


= 2{k l k 2 ) 








h 


= 2{k lPl ) 


i'i 


= 2(k 2Pl ) 


s 4 


t 2 


= 2{k 2 p 2 ) 




= 2(k lP2 ) 


•53 


t' 


= 2(k 2 k 3 ) 


v! 


= 2(hk 3 ) 


«5 



l(pm) +m( + mi (B.9) 
All momenta are chosen incoming: k\ + k 2 + k 3 + p\ + p 2 = 0. Only five invariants are independent: 

s 3 = - s 5 - t 2 - u 2 - m\ + m\ 

54 = — S5 — t\ — u\ — m\ + m\ 

55 = + S + t' + v! 

u = — s — t\ — u 2 

t' = - s-t 2 -m (B.10) 

From the relation 54 = s + t 2 + u 2 — m\ + m\ one can see that the limit S3, S4 — > is equivalent to 
the Born kinematics. The differential cross section in the Born approximation therefore includes a factor 



o d 2 a B Tr tv (An) 2+e f t 2 u 2 - sml\ e r . , \-^,,, R |2 ^ „ 



'As usual A(a;, y, z) = x 2 + y 2 + z 2 — 2(xy + xz + yz) 

"Note that in the formulae (D. 11) and (D. 12) of [p5|] the indices of the mo 1 term have to be interchanged 
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with the n dimensional spin-color averaging factors for different strongly interacting incoming states 

K ^ = W^ K " = 4(1 - e) 2 (iv 2 _ 1)2 K i9 = K m = 4(1 _ e)iV(iV 2 _ i) < B - 



The matrix element for the real gluon emission corresponds to a differential cross section in four vari- 
ables. 



s' 



dMil " 2r(l - 2 £ ) V ^ 2 J (-4 + m?) 1 - /dn> ' |A< 1 (B " 13) 



The total partonic cross section is defined as 

a= / — — = s + t 2 + m 2 2 -m\ + — ± 

Jtf in JO at 2 dSi t 2 

min/max 



/ _ s + ml-m\± A 1 / 2 (s, m\,m 2 2 ) ^ ^ 



To integrate over the angle of the final state gluon analytically, one chooses different parameterizations 
of the phase space. The angular integration is performed in the center-of-mass frame of p\ and k 3 . One 
of the three dimensional components of k 2 , k\ or p 2 is taken parallel to the z axis. In the first case the 
momenta are [ 38 ] : 

ki = (-wi, 0, -ps^, -pc^ + io 2 ) 

k 2 = (-w 2 ,...,0,Q,-w 2 ) 

h= (w 3 , ...,w 3 sis 2 ,w 3 sic 2 ,w 3 ci) 

Pl = (E 1 ,...,-W 3 SlS 2 ,-W 3 S 1 C 2 ,-W3Cl) 

p 2 = (E 2 , ...,0,psip,pcip) (B.15) 

Si, Ci are defined as the sine and cosine of the angles 9i,0 2 ,i^. The angles 9j are connected to the angular 
integration for the additional gluon: 

r /"it 

n-3 / „n-4 



dn = / d9 lS ^- A / dQ 2 s n 2 ~* (B.16) 
Jo Jo 

The non-invariant variables can be expressed in the invariants: 

s + u 2 s + t 2 s 4 

wi = — , w 2 = — . w 3 



2\J S4 + m 2 2a/s4 + m 2 %\J~S4 + m\ 

s 4 + 2m? £ 2 + u 2 + 2m 2 , 
— — , ,> E 2 — - 

y/ (t 2 + U2) 2 - 4sm| t 2 (s4 - m| + mf) - s(u 2 + 2m 2 .) 

P = — C0 = (d . 1 / j 

2y / s A + m\ (s + t 2 )^J{t 2 + u 2 ) 2 -Asm 2 2 

Inserting these equalities into the matrix element after partial fractioning leads to the typical integral [ [38| ] 

dV- tt k,leN (B.18) 

{a + bdf (A + B Cl + Cs lC2 ) 1 

Note that it is not necessary to perform this integration analytically, if the divergent regions of phase 
space are regularized e.g. by the subtraction method [Eq]. 
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Cut-off Method 



The radiation of on-shell gluons leads to divergences in the phase space integration in two limiting cases: 
(i) the gluon is soft, all components of its four vector vanish k^ — > [fx = 0, 1, 2, 3]; (ii) the gluon is 
collinear to another massless particle with the momentum p, i.e. k ~ p + k± with k± — > 0. In both cases 
the invariant vanishes, (pk) — ► 0, leading to infrared divergences in the matrix element. In the limit of a 
soft gluon the invariants for the three body process approach 



h + m 1 



s 3 -> f -» i 2 

S4 — ► u' — » U2 — > Ui + mf — m\ s§ — > s (B.19) 
The integration of real gluon matrix element is split into two regions, corresponding to soft and hard 



gluons [38]: 



ds 

4 C?t2^S4 



+ 



cIsa 



dsA 



A/ 

d 2 a 



d 2 a 

dt2dS4 



dt2dS4 



approx 



ds ^ 
& 4 dt2dsA t 



(B.20) 



In the second integrand the limit A — > has to be checked numerically. The soft gluon matrix element 
in the first integrand is evaluated in the eikonal approximation, where the gluon momentum is neglected 
compared to any other variable [<J(s4)]. The angular integral is analytically evaluated using the angular 
integrals eq.( B.18 ). In addition to soft poles in e logarithms log J A[j = 1,2] appear. These have to 
be added to the final expression. To cancel the A dependence of the hard gluon part, these terms are 
rewritten as 



ds4, log ( 



io g (sr/// 2 ) 



A 



1 

S4 



6(34) = J dsa. 

In addition to the soft singularities in the gluon emission matrix element collinear poles arise in next-to- 
leading order. As long as the final state particles are massive, these have to be absorbed completely into 



(B.21) 



the renormalization of the partem densities, as described in chapter 2.2.2. The additional terms are of the 
form 



d 2 a hIF 



dt2ds4 2ir 



IE + log 



/ 47T/X 2 



X -Pu{x) is 2 

x V 



dt2dS4 



(B.22) 



xki 



The splitting function P{j describes an incoming parton i changing to j^] 

1 



Pgg(x) = 2C A 

C F 



X 



Pgq(x 



1 

1 h (1 - 



+ 



X 



l + x(l 



+ 5(1 



11 



Ca - ^T R nf 



Pgg(x) = T R (X 2 + (1 - Xf) 



P qq (x) = C F 



1 + x 2 
1 — X 



(B.23) 



3 The color factors for the SU(3) are N = 3, Ca = 3, Cf = 4/3, T R = l/2 
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with the + distribution defined as 

(F(x)) + = F(x) e(l-x-P)-5(l-x- p) C P diF{i) (B.24) 

J o 

in the limit (3 — ► 0; the parameter separates soft from hard gluons. Numerically the + distributions can 
be evaluated following 



J dx [f{x) (L(x)) + + g(x)) = J dx lf(x)L(x) + g(x)) - /(l) 



1-/3 r l-l3 

dx 



L(x)+ 1 



l-T 

(B.25) 



where L = J Q T dxL{x) has to be calculated analytically and (3 — ► 0. The second part of the + distribution 



eq.( |B.24| ) contributes to the soft gluon part of the splitting function. After performing the integration 
over the momentum shift x using the 5(1 — x) term, the result follows the Born kinematics and cancels 
the collinear divergences of the virtual corrections. The separating parameter can be linked to the cut-off 
A via (3 = A/(s + £2) and /? = A/(s + U2) for a shifted momentum &2 or k±. 

For hard gluons the four momentum conservation 5(s4) included in the Born differential cross section 
can be used to remove the integration over the shift in the four momentum, i.e. the mass factorization 
term cancels the collinear divergences arising from the virtual and from the real corrections. 



Subtraction Method 



The subtraction method [56] is used to calculate the neutralino production cross section, i.e. there are no 
massive strongly interacting particles present in the final state and qq is the only incoming state. 

Given a divergent real gluon matrix element, a subtraction matrix element a A,s is constructed for the 
three particle final state, in order to remove the soft and collinear singularities point-wise from the gluon 
emission phase space. 



d 2 a R 2 d 2 a A > 3 
dt2ds4 dt2ds4 



„2 g °_ s 2 d a " . . S^\\stR\ 2 T.13,2 -r.23,1 



J2\ MR \ -P 13 ' 2 -^ 23 ' 1 (B.26) 



In case of only two initial state partons k\ , hi and one emitted gluon ^3 the dipole terms 

have to be added to ^ |A^ R | 2 . The rescaling factor for the split incoming momentum is given as 
{k\k,2 + k^k\ + k^k2)/k\k2 and re-defines the kinematics of the Born matrix element k'^-p'- e.g. for V 

k-^ — x k\ 
k' 2 = k 2 

, 2p j K + 2p j K' , 2piK 

P'i = Vi ~ [ K + K,y (K + K^ + ^K (B.28) 



x = 

13,2 



68 



APPENDIX B. RADIATIVE CORRECTIONS 



with K = k\ + k 2 + &3 and K' = k[ + A;^. The result for V 23 ' 1 is obtained by interchanging the incoming 
quarks. The kernels V 13,2 and V 23,1 are of the form 



V q9 ' q {x) = 8?r/x 2e a s P qq {x) = 8TT^ e a s C F 



2e. 



X 



-(l + x) + e(l-x) 



(B.29) 



Since this subtracted real gluon emission matrix element is finite, da R and da A ' 3 can be evaluated in 
n = 4 dimensions. However, this does not hold for the virtual corrections matrix elements, where one 
has to use the dipole terms da A ' 2 , which arise from the exact integration of the subtraction term da A ' 3 
over the additional gluon emission phase space. This results in soft and collinear poles in e, therefore the 
dipole moments for the calculation of the virtual subtraction term have to be evaluated in n dimensions. 
The integrated soft and collinear gluon subtraction yields 



dU 



dx I(x, e) 



,da B 
~dh 



+ 8 



xk\ 



,da 



B 



dt 2 



xk2 



(B.30) 



The integration of V q9 > q (x) and the mass factorization leads to the integration kernels for quark-antiquark 
and quark-gluon incoming states: 



Pis 

2^ 



C F ( 4TT/I 



T(l-e 

xs 



,2\ e 



5(1 -x) 



-P qq ( x )\og[®F\ +C F [4 



1 3 
log(l - x) 



X 



fqgix) = 



a s 
2vr 



- P Q9 log 



Q] 



(1 — x) 2 s 



+ 2T R x(l-x) 



- 2(1 - x)log(l -x) + l-x 

(B.31) 



Since using the subtraction method the gluon emission angles are integrated numerically the mass fac- 
torization terms are not added analytically to the real gluon matrix element. Instead of using the S(s4) 
distribution to remove the integration over the momentum shift x, it is kept for the phase space integra- 
tion. The convolution is then performed numerically. The separation into virtual-soft and hard gluon 
contributions is by no means unique. 



B.2. Hadronic and Differential Cross Sections 

For a general hadronic cross section the partons are treated as parts of a massive hadron hj[j = 1, 2]. 
The hadronic variables are referred to as capital letters. 

hi(K!) + h 2 {K 2 ) - Xi(pi) + X 2 ( P2 )[+g(k 3 )} (B.32) 

The additional gluon is again assumed massless, the pj [j = 1, 2] will in general have different masses 
mi,m 2 . We introduce the invariants: 



S = 2(K 1 K 2 ) T x = 2(K lPl ) XJ X = 2(k 2Pl ) 

T 2 = 2{K 2P2 ) U 2 = 2(K lP2 ) (B.33) 
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All momenta are chosen incoming. After integration of the gluon emission angles the real emission 
matrix element together with the parton densities is a differential cross section with respect to four 
variables including some general function F: 

ctot = / dx\ / dx 2 / dt 2 / ds± F 



o 



min _ ( m l + m 2) 
Xi — - 



-mm 



AS 

(mi + TO2) 2 



45xi 

sf** = s + t 2 + mi,-m\ + 



,2 



,min/max s + m^ - mf ± A 1 / 2 (s , mf, 
J 2 — ^ (o.34) 

This can be rewritten for the on-shell subtraction in the 54 channel, when the final state particle X\ is 
produced via decay of one on-shell squark [||] : 

rl rl rsf** rtf** 

ftot = / dxi / dx 2 / ds4 / F 







s + m 2 — m^ — 2-v / sm^ 



^min/max s - s 4 + m| - m\ ± y^s - s 4 + m 2 , - mf ) 2 - 4sm| 
t 2 = 

The integration variables x\,x 2 are only used to compute the total hadronic cross section. Therefore the 
subtraction in S3 involving the particle X 2 can be obtained by exchanging p\ <-> p 2 in the subtraction 
matrix element and in the phase space. This holds only if the gluon emission angle d£l is integrated out 
completely. 

To be able to compute differential cross sections with respect to the transverse momentum of one of the 
final state particles or the rapidity 

2 T 2 U 2 - Sml t 2 u 2 - sm\ 1 / T 2 \ 

p t = - s = y =2 log \Th) (B36) 



the integration has to be ordered differently 

„ e max „ 1 

X\S + T 2 



o-tot = / dp T I dy I ds A I dxi 2p T S F 



j mlx JO Jx 



max = J_ A 1/2 (S 2 2) 



.max 



arcosh 



S + m 2 — m\ 
2y/s(p*+ml) 
s ™*= S + T 2 + U 2 + m 2 2 -m\ 
. S4 _r 2 -m 2 + m? 
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with x 2 = (54 - x\XJ 2 - ml + mf) / (x±S + T 2 ). 

The subtraction of differential cross sections in the S3 channel requires a substitution of one of the two 



angular integrations in dft, which is implicitly included in eq.( |B.37| ). We refer to the integration 
boundaries therein clS S ^ . X ^ . 



'4 fl 
At 



x* JO JO 



d01 



d6 2 s 



2*1 



- s max j .jjiai 

/ dS3 / cfei / C?S4 / 

JO A™ in is™° JO 



2(s 4 + mj) 
S4 y s — S4 — mf + 777,3 



(B.38) 



The integration borders are 



,r-i 



7777 



2(5 + T 2 + ?7 2 + ml) 
1 

! 2(5 + U 2 ) (m 2 2 S + m2?7 2 + s 3 U 2 ) 

.22, 



T 2 -U 2 - 2m 2 2 + \A T 2 + ^) 2 - 4777IS 



2m 2 S + 2m\m\S — 2m 2 Ss 3 — Ss 3 — 2m 2 ST 2 



^min/max 



Ss 3 T 2 - 2m 2 U 2 + m 2 m{U 2 - 3m 2 s 3 U 2 + m{s 3 U 2 - s 3 U 2 - 2m 2 T 2 U 2 - 2s 3 T 2 U 2 
s 3 (±m 4 2 S 2 - Am 2 2 m{S 2 + Am 2 2 S 2 s 3 + S 2 sl + Am 2 2 S 2 T 2 + 2S 2 s 3 T 2 + 5 2 T 2 2 + 4m\SXJ 2 
- 4mlm\SU 2 + Qm 2 2 Ss 3 U 2 - 2m\Ss 3 U 2 + sSs 2 3 U 2 + 2m 2 2 ST 2 U 2 + 2m\ST 2 U 2 

+ 2Ss 3 T 2 U 2 + m\ul - 2m 2 2 m 2 Ul + mfU^ + 2m 2 2 s 3 Ul - 2m\s 3 U 2 2 + s 2 3 U^ ' " 



•5.3 



777 2 T 2 — m\T 2 — s 3 T 2 — 2m\x\S — x\Ss 3 — x\SU 2 



' 2{ml Xl S + m 2 2 T 2 + s 3 T 2 ) 

=f ^(-m 2 T 2 - m\T 2 - s 3 T 2 - 2?77, 2 xi5 - x\Ss 3 - x\SU 2 ) 2 

1/2 



- 4m 1 (xiS + T 2 ){m 2 xiS + m 2 T 2 + s 3 T 2 ] 



xi{S + U 2 ) + T 2 + m 2 



777,1 



for s 3 < S3 
for S3 > S3 



,* _x\S + T 2 + x\U 2 + m 2 — m\ 
' 3 2(xi5 + T 2 + xi[/ 2 + m2) 



T 2 - Xl U 2 - 2m 2 -J(T 2 + Xl U 2 f - 4m\x x S 



(B.39) 



B.3. Scalar Integrals 
Dimensional Regularization 

The analytical expressions for the virtual corrections contain scalar integrals which are multiplied by 
polynomials including the Mandelstam variables. The scalar integrals A, B,C, D in case of dimensional 



regularization Q35| , p7Q are defined in n = 4 — 2e dimensions and are used to calculate the production 
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processes of neutralinos/charginos and stops. 



C{{ Pl };{m 3 }) 



4— n 



i(27r) n [q 2 - m\] [(q + P1) 2 - m\] [(q + p 12 ) 2 - m§] 



d n q 



D({ Pl };{ mj })-- 

Pijk =Pi + Pj +Pk 



4— n 



i(2vr) n [q 2 - mf] [(q + pi) 2 - m|] [(g + pi 2 ) 2 - m|] [(g + P123) 2 - m 2 ] 



(B.40) 



The definition of the one and two point functions A, B follows the same conventions. Using this in- 
tegration measure, the non-absorptive virtual contributions are real and the integrals have got integer 
dimension. The infrared and collinear divergences occur as poles l/e k [k = 1,2] [M g). p is the 
renormalization scale of the process. 

As long as only CP conserving observables are calculated, the typical combinations of couplings in front 
of the scalar integrals are real, therefore only the real part of the scalar integrals is needed. If, as described 
in appendix |aJ, we chose the parameters in the Lagrangean as being complex we also need the imaginary 
parts of the scalar integrals. 



The expressions for finite integrals have been taken from the literature [58, 55]. The divergent integrals 



have been calculated using Cutkosky cut rules and dispersion relations [59], and most of them are also 
present in the literature [§]. 

A typical singular scalar three point function occuring in the neutralino/chargino production cross section 
is: 



Cipukr, m x ,0,0) 



t x = t 
M 2 



rn , 



171; 



Ce 

h 



log 



+ log 2 
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+ Li 2 
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Li 2 
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mi 
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+ log ( 5 ) log 



-t x 

M 2 



(B.41) 



The definition of the momenta and Mandelstam variables follows appendix |B.1[ The factor in front 
contains the typical MS terms 



a 



16vr 2 



/ 47T^ 



2\ e 



(B.42) 



The scale M is either chosen as the mass of the outgoing particle or, for different masses in the final 
state, as the averaged mass. The factor C e is a common factor of all virtual gluon contributions and can 
be pulled out of all renormalization contributions, i.e. it occurs as an over-all factor in the regularized 
and renormalized matrix element and can then be evaluated in the limit n = 4. 



72 



APPENDIX B. RADIATIVE CORRECTIONS 



One scalar four point function is 



D(k 2 ,k 1} px; 0,0,0, m x ) 
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(B.43) 



The log(l + M^M| / (sfRj)) multiplied with the terms in brackets originates from the analytical contin- 
uation of the Dilogarithm ||59|]. This method of calculating scalar integrals omits the typical roots, which 
appear by using the Feynman parameterization after partial fractioning. Thereby, large cancelations are 
absent in our results, which improves the numerical accuracy. 



Massive gluon regularization 

For the massive gluon regularization scheme, as it is used for the calculation of the decay processes, the 
conventions concerning the scalar integrals are exactly the same as for dimensional regularization in the 
limit n = 4. In this case we obtain divergences in form of logarithms of the gluon mass log A 2 [59]. 
For linear divergences the log A 2 description can naively be translated into 1/e, whereas for higher diver- 
gences this is more involved. 
For the gluino decay the divergent scalar integral 



C(pi,P2', mi, A, m 2 ) 



C e Xe 



mim 2 (l — x 2 ) 



log(x s ) 
+ Li 2 

+ 2 log 2 



log 



(— ) 

\mim 2 J 



^log(x s ) + 21og(l-x 2 ) 



x s — + Li 2 
m 2 



mi 
m 2 



)+Li2 (l-*^)+Li 2 (x 2 ) 



(B.44) 



is used in the massive gluon regularization scheme. The factor C e is taken in the limit n = 4, and 



Pi 



J- s 



ml 



Am\m2 
s— (mi — r«2) 2 



m 2 
- 1 



4mim-2 
s— (mi — 1T12)' 1 



+ 1 



(B.45) 
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For the stop production the same integral is needed in the case of equal masses and regularized dimen- 
sionally. There the divergent logarithm has to be replaced by a linear pole log A 2 — ► 1/e + log ji 2 R 



C(pi,p 2 ;m 1 ,0,m 2 ) 



m\mi(\ — xf) 









\og{x s ) 


-- - log 


( M l ) 




e 


\m 1 m2 J 



-log(x s ) + 21og(l-x 2 ) 



(B.46) 



The same relation holds for the other scalar integrals used calculating stop decay widths. 
B.4. Counter Terms 

Renormalization of the external masses in the on-shell scheme preserves the Slavnov-Taylor identity [ [38| ] 

kl M^ v = Mf host oc k 2v (B.47) 

where M.^ is the matrix element for the production of two massive SU(3) charged particles in gluon 
fusion and Aif host the ghost contribution, both of which are present in LO and NLO. This identity has 



been used to calculate the ghost contribution of the stop pair production, as described in chapter |1.2.4 . 
The mass counter terms for a heavy quark t and a squark q are 
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The strong coupling is renormalized in the MS scheme 
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-JV 



— log — |^ - - log -4- 
12 % 2 3 



-n 



(B.49) 



is the renormalization scale of the process, rif = 6 the number of quark and squark flavors. The 
j3 function gets contributions from the gluons, the gluinos, the quarks, and the squarks. For the MSSM 
the sign of this coefficient is the same as for the Standard Model. To decouple the heavy particles from 
the running of a s the massive log to 2 terms have to be subtracted i.e. the strong coupling constant is 
effectively evaluated as the usual low energy Standard Model QCD coupling. 
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In addition to these Standard Model like renormalization constants, which remove the UV divergences 
from the virtual correction matrix element, finite counter terms have to be added to restore the supersym- 
metric Ward identity, as described in chapter 1.5. These are no counter terms observables arising from 
divergent vacuum fluctuations in field theory, i.e. even the weak coupling will contain a new counter term 
oc a s although physically Gp should stay unrenormalized in supersymmetric QCD. 



C. Analytical Results for the Stop Decay Width 



As an example, the analytical results for the decay of a light stop t\ into top and gluino are given. The 
decay width in next-to-leading order may be split into the following components: 

r NLO = r LO + Re (Ar t + Ar § + Ar n + at,, + Ar r + Ar c + Ar, + Ar dec ) (c. i) 

To allow for more compact expressions we first define a few short-hand notations: 

A 1 / 2 (m? , m~,m 2 ) 

Vabc = m a + m b -m c a 2 § = m t m- g sin(20) N = 167r ^3 

ti 



where a,b,c = g,t,j with j representing tj. The different contributions defined in eq.( |C.l[ ) are listed 
below: 



lowest-order decay widthQ 

r LO = 8NC F 7ra s [-fjt m + 2a 2§ ) N = N\M B \ 2 (C.3) 
top self-energy contribution: 

2NC F -Ka s \M B \ 



AT 4 = ™ ^f*ocs\JVIb\ f 2(1 _ e)A{mt) + 2A{ } _ A{ } _ A{ j 
ml y 

+ ii t r g B(p t ;mg,m ii ) + fi t2 gB(p t ; m g , mQ 

~ 4m i a 2§ [ B (Pu m g, ™>ii ) ~ B (Pf, m g, m t 2 )] 

+ 2mf [ng t iB(p t ;m~g,m~ tl ) + fj,g t 2B(p t ; m g , mjj - 4m 2 B(p t ; A, m*)] | 

^i{-4(m t - 2 ) -^4(m fl ) - ^ gt iB{p t ;m g) m- tl ) + n gt2 B(p t ; m g , m ?2 )| 

(C.4) 



16AAAC|7r 2 a^ 2 



where the definitions of the scalar integrals are given in eq.( |B.40| ). When renormalizing external masses 
in the on-shell scheme the scalar function B(p; m a , m^) = dB(p; m a , m^/dp 2 appears. 



'The Casimir invariant for the gauge group SU(3) is Cf = 4/3 
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gluino self-energy contribution (for nj = 6 quark flavors): 

Ar s = ( n f- 1 ){- A ( m q) + { m \ + m$)B(pftmfaQi) + 2m 2 g {m 2 g - m~)-B(p§; mj, 0)| 

+ ^y^i^l^J |2A(m t ) - A(m fl ) - A(m ?2 ) + ^^(p^m^mt) + ii 2 gtB(pg;mi 2 ,m t ) 
- 4m~a 2§ [B(pg\ , m fl , m t ) - B(p g ;m i2 , m t )] 
+ 2m 2 g [fi gtl B(p g ;m il ,m t ) + li g nB{p g \ m t ~ 2 , m t )] j 

+ ^^l^| 2 { (1 _ 6)A(rn .) _ 4m jB(p g ; A ,m 5 )} (C.5) 

diagonal stop self-energy: 

ATu =8C , FAA7ra s |A4 jB | 2 |5(p i - i ;mg,mj) - B(p ii ;X, mjj + 2a 2 gB(p ii ;m g ,m t ) 

- ngtiBiPh ; m 9> m *) - 2m | ^(p^ ; A > m tx ) } (C- 6 ) 

[The off-diagonal mixing contribution 



128AAA^C|vr 2 a2 2 J 4m 2 m| 
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AFl2 = ^2 _i 2 S *2* c 2e A K~J " A ( m ti) + — — B (Pti ! ™s> mt) (C.7) 



rrt~ — m~ ^" <T n 



is absorbed into the renormalization of the mixing angle for 6 = 6(mf ), as described in chapter 1.2.4 ] 
vertex corrections: 



AT V = 64MTT 2 a 2 s NCl F[ + a 2§ Ff + a\^l 



+ 32Nir 2 a 2 s N 2 C F Ff + a^F? + a 2 § F^ + 8Nira s \M B \ F B (C.8) 
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Fi = - 2en~gtiB(pt i ;m~ g ,m t ) + 4(m 2 - m|)JB(p § ; A, m§) 

+2m 2 [B(p t ;m g ,m i2 ) - B(p t ;rrig,m il )] + Am 2 g B(p g ;m t ,m h ) 
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corrections from real-gluon radiation: 
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renormalization of the coupling constant as defined in eq.( p.49| ) 
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finite shift of the Yukawa coupling relative to the gauge coupling in MS, as described in chapter |L5 
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decoupling of the heavy flavors from the running strong coupling constant: 
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